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ABSTRACT 


A  numerical  study  is  performed  to  test  the  hypothesis 
that  the  Pacific  polar  low  forms  primarily  through  baro- 
clinic  instability./  Linear  growth  rates  computed  for  a 
jet  profile  typical  of  conditions  during  the  formation  of 
Pacific  polar  lows  yields  a  growth  rate  maximum  at  1400  km. 
Although  the  horizontal  scale  of  this  disturbance  agrees 
with  observations  of  the  Pacific  polar  low,  this  disturbance 
is  shallower  and  forms  further  north  relative  to  the  jet 
than  indicated  by  observations.  Nonlinear  simulations  of 
the  1400  km  disturbance  indicate  that  the  effects  of  convec¬ 
tion  and  surface  heat  and  moisture  fluxes  due  to  initially 
large  air-sea  temperature  differences  are  needed  for  proper 
vertical  structure,  storm  position  and  intensification. 
Energy  budgets  of  model  output  indicate  that  the  wave  grows 
through  baroclinic  instability  and  that  the  diabatic  effects 

enhance  the  baroclinic  conversion  within  the  storm. 
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The  second  jet,  referred  to  as  the  genesis  jet,  is 
typical  of  conditions  during  the  genesis  period  of  north¬ 
eastern  Pacific  polar  lows  (Mullen,  1979).  Mullen  (1979) 
defined  this  period  to  be  "at  the  latest  hour,  either  0000 
GMT  or  1200  GMT,  before  the  cloud  mass  in  question  showed 
any  discernable  comma  or  hook  shape  in  the  infrared  (IR) 
imagery."  Mullen  composited  NMC  analyses  for  22  north¬ 
eastern  Pacific  polar  lows.  In  all  cases,  the  width  of 
the  cloud  mass  at  the  genesis  period  was  at  least  500  km. 

It  was  expected  that  this  jet  would  be  appropriate  for 
studying  the  polar  low. 

The  mean  winter  and  genesis  jets  of  Lau  (1973)  and 
Mullen  (1979)  were  transformed  to  model  coordinates  in  the 
following  manner.  First,  the  channel  width  of  the  model  was 
chosen  to  be  6400  km  with  the  jet  core  in  the  center  of  this 
channel.  Second,  at  all  mandatory  levels  up  to  300  mb  in 
the  genesis  jet  and  up  to  100  mb  in  the  mean  winter  jet, 
wind  speeds  were  tabulated  in  250  km  increments  over  an 
interior  2500  km  meridional  slice.  It  was  assumed  that  the 
zonal  wind  for  both  jets  decreased  linearly  to  zero  at  both 
the  north  and  south  walls.  The  shape  of  the  genesis  jet 
above  300  mb  was  assumed  to  resemble  tnat  of  the  mean  winter 
jet.  Third,  these  wind  speeds  were  linearly  interpolated 
on  each  mandatory  pressure  surface  to  the  model  horizontal 
grid  points.  Fourth,  one  pass  with  an  81  point  smoother- 
desmoother  (Shapiro,  1  9  7  0  )  was  applied  to  these  fields  to 


25 


II.  NUMERICAL  PROCEDURES 


The  NORAPS  forecast  model  has  several  features  which 
make  it  applicable  for  extensive  numerical  simulations. 
First,  the  domain  size  and  the  horizontal  and  vertical 
resolution  in  the  model  may  be  varied.  Second,  dynamics  and 
physics  appropriate  to  subsynoptic  features  are  included  as 
described  in  Appendix  A.  Exact  values  for  various  param¬ 
eters  are  given  in  the  description  of  each  model  experiment. 
A  map  factor  of  1  is  used  over  the  entire  domain  for  all 
forecasts.  Also  the  Coriolis  force,  f,  is  calculated  using 
the  6-plane  approximation, 

f  =  fo  +  s0(y-yQ)  (2.1) 

A.  INITIAL  DATA 

1 .  Initial  Jet  Prof i 1 es 

Two  baroclinically  unstable  jet  profiles  were  chosen 
for  this  study.  The  first  jet  describes  the  mean  winter 
conditions  over  the  Pacific.  The  conditions  were  derived  by 
Lau  (1978),  who  composited  National  Meteorological  Center 
( N  M  C )  analyses  from  850  to  100  mb  for  11  winters  in  the 
sector  from  150E  to  150W  and  from  20N  to  70N.  Mullen  (1979) 
included  1000  mb  winds  estimated  by  Crutcher  and  Meserve 
(1970).  It  was  expected  that  this  jet  (referred  to  here¬ 
after  as  the  mean  winter  jet)  would  be  appropriate  for 
studying  typical  oceanic  sy no p t i c - sc  a  1 e  lows. 
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also  shows  how  barotropic  and  baroclinic  instabilities 
contribute  to  the  generation  and  destruction  of  the  eddy 
kinetic  energy  of  the  most  unstable  waves. 

The  second  major  section,  Chapter  IV,  describes 
nonlinear  model  simulations  of  the  polar  low.  Detailed 
analyses  of  wave  structure  and  energy  budgets  are  presented 
for  adiabatic  and  diabatic  simulations.  The  sensitivity  to 
the  air-sea  temperature  difference  is  discussed. 
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To  test  this  hypothesis,  the  numerical  model  used  in  the 
Fleet  Numerical  Weather  Center's  Navy  Operational  Regional 
Atmospheric  Prediction  System  (NORAPS)  is  used  (Appendix  A). 
This  model  includes  convective  and  non-convective  precipita¬ 
tion  processes  and  state-of-the-art  planetary  boundary  layer 
(PBL)  physics.  In  addition,  the  model  may  be  run  with 
various  combinations  of  grid  spacing  and  resolution.  This 
makes  the  model  quite  useful  for  studying  waves  of  various 
sizes. 

Interpretation  of  model  results  is  accomplished  through 
budgets  of  kinetic  and  available  potential  energy  as  given 
in  Robertson  and  Smith  (1983).  These  equations  (Appendix  B) 
are  split  into  zonal  and  eddy  components  and  the  exact  form 
of  the  available  potential  energy  is  used.  Since  the  fore¬ 
cast  model  uses  a  closed  domain,  the  eddy  components  can  be 
directly  related  to  a  single  wave  in  the  forecast  domain. 

This  study  consists  of  two  major  sections.  The  first, 
Chapter  III,  describes  linear  growth  rates  computed  using  a 
linearized  version  of  the  NORAPS  model.  Two  sets  of  initial 
conditions  are  used.  The  first  is  a  mean  winter  jet  and  the 
second  is  a  jet  typical  of  conditions  during  the  development 
of  Pacific  polar  lows.  Small  perturbat i ons  are  introduced 
into  the  initial  fields  and  an  adiabatic  version  of  the 
model  is  integrated  while  restoring  the  initial  mean  zonal 
fields  every  time  step.  The  effects  of  linearized  friction 
and  wave-zonal  flow  interaction  are  discussed.  This  section 
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Gall  (1976c)  demonstrated  that  convection  can  serve  to 
increase  the  vertical  extent  of  medium-scale  (2000  km) 
waves . 

Surface  fluxes  of  heat  and  moisture  also  must  be 
considered  since  large  air-sea  temperature  differences  are 
often  associated  with  Pacific  polar  lows  (Mullen,  1979). 
Large  upward  fluxes  of  heat  can  reduce  the  low-level  static 
stability  which,  as  explained  above,  can  result  in  baro- 
clinically  unstable  1000  km  waves.  The  role  of  surface 
fluxes  in  strengthening  disturbances  has  been  discussed  by 
Petterssen  (1956),  Danard  and  Rao  (1972)  and  Forbes  and 
Merri  tt  (1984). 

Based  on  the  observational  and  analytic  studies 
described  above,  the  major  hypothesis  of  this  study  is  that 
Pacific  pojar  lows  are  formed  primarily  through  baroclinic 
instability.  They  are  a  preferred  mode,  especially  in 
winter,  because  of  the  unique  vertical  temperature  profile 
forced  by  a  polar  air  mass  moving  over  relatively  warm 
water.  The  upward  sensible  heat  flux  in  the  polar  air 
reduces  the  low-level  static  stability  and  therefore,  R.. . 
This  is  precisely  what  is  needed  to  set  up  barocl i ni cal  1 y 
unstable  waves  on  the  scale  of  the  Pacific  polar  low 
according  to  the  theoretical  studies  discussed  above.  It 
is  also  hypothesized  that  surface  fluxes  and  latent  heat 
release  through  convection  modify  the  baroclinic  wave  by 
further  intensifying  it  and  increasing  its  vertical  extent. 


tends  to  shift  the  maximum  growth  rate  to  shorter  wave¬ 
lengths  (Gambo,  1970a;  Tokioka,  1970;  Mansfield,  1974;  Gall, 
1976a;  Duncan,  1977;  Staley  and  Gall,  1977).  However,  the 
role  of  baroclinic  instability  is  still  in  doubt  since  short 
baroclinic  waves  tend  to  be  quite  shallow,  as  opposed  to  the 
observed  deep  structure  of  the  Pacific  polar  low.  In  addi¬ 
tion,  Sardie  and  Warner  (1983)  used  a  quasi -geostrophi c 
model  to  show  that  the  maximum  growth  rate  for  conditions 
typical  of  those  during  the  formation  of  the  Pacific  polar 
low  is  found  at  a  wavelength  of  2700  km. 

The  presence  of  convection  suggests  that  it  may  be  a 
significant  factor  in  the  development  of  the  Pacific  polar 
low.  Rasmussen  (1979)  has  shown  that  CISK  is  capable  of 
producing  "m i n i - hu r r i c anes"  similar  to  the  Atlantic  polar 
low.  However,  CISK  requires  a  pre-existing  circulation 
which  must  be  explained.  Also,  CISK  requires  a  small 
vertical  wind  shear,  a  condition  not  often  met  during  the 
development  of  the  Pacific  polar  low  (Mullen,  1979). 

Therefore,  it  appears  that  some  combination  of 
baroclinic  instability  and  diabatic  processes,  such  as 
convection,  are  responsible  for  the  development  of  the 
Pacific  polar  low.  It  has  been  shown  that  Rj  coupled  with  a 
CISK-type  latent  heating  shifts  the  most  unstable  wavelength 
from  the  synoptic-scale  to  1000-2000  km  waves  (Gambo,  1976; 
Rao  and  Ferreari,  1979;  Sardie  and  Warner,  1983).  Also, 


surface  (Reed,  1979).  The  Pacific  polar  low  tends  to  form 
north  of  the  polar  jet  in  a  strong  baroclinic  zone  within 
a  polar  air  mass.  The  overall  appearance  of  this  system 
resembles  that  of  typical  synopt i c- seal e  extratropi cal 
systems,  but  of  smaller  scale. 

The  purpose  of  this  study  is  to  investigate  the 
reasons  for  the  growth  of  the  Pacific  polar  low.  It  has 
been  proposed  by  Reed  (1979)  and  Mullen  (1979)  that  either 
barotropic  instability,  Conditional  Instability  of  the 
Second  Kind  (CISK)  or  baroclinic  instability  may  be  respon¬ 
sible,  either  singly  or  in  some  combination.  Although  the 
conditions  for  barotropic  instability  are  often  met  in  the 
Pacific  during  polar  low  formation  (Mullen,  1979),  Gall 
(1976a)  showed  that  barotropic  effects  are  about  an  order  of 
magnitude  smaller  than  baroclinic  effects  for  wavenumber  15. 

Many  analytic  studies  support  the  argument  that 
baroclinic  instability  is  responsible  for  the  growth  of 
the  Pacific  polar  low.  The  earliest  studies  of  baroclinic 
instability  (Charney,  1947;  Eady,  1949)  indicate  that 
4000  km  waves  are  most  unstable  and  observations  verify  that 
this  scale  dominates  in  extratropi  cal  regions.  However, 
more  recent  analytic  studies  show  that  the  growth  rate 
spectrum  is  very  sensitive  to  the  Richardson  number  (R-j), 
which  is  the  ratio  of  the  static  stability  to  the  square  of 
the  vertical  wind  shear.  Reduced  R:  in  the  lowest  levels 


I.  INTRODUCTION 


The  availability  of  global  satellite  photographs  on  a 
daily  basis  has  exposed  meteorologists  to  weather  disturb¬ 
ances  which  were  previously  not  known  to  exist.  In 
particular,  systems  called  "polar  lows"  (Harley,  1960) 
have  been  found  to  exist  primarily  over  the  oceans  in  the 
winter  seasons  (Reed,  1979;  Rasmussen,  1979).  These  storms 
which  have  horizontal  scales  of  only  500-1500  km,  can  be 
quite  violent  and  occasionally  exhibit  thunder,  lightning, 
hail  and  hurri cane-force  wind  gusts  (Stevenson,  1968; 
Suttie,  1970;  Lyall,  1972). 

It  has  been  suggested  that  two  distinctly  different 
types  of  polar  lows  exist  (Rasmussen,  1981;  Locatelli 
et  a  1  . ,  1982;  Sardie  and  Warner,  1983).  The  first  of 
these  is  the  Atlantic  polar  low,  which  typically  has  a 
spiral  structure  and  occasionally  an  "eye"  similar  to  that 
found  in  tropical  storms.  This  disturbance  exhibits  a 
horizontal  scale  of  500-1000  km  and  maximum  amplitude  at 
the  surface.  In  addition,  it  tends  to  form  within  a  polar 
air  mass  north  of  the  polar  jet  in  an  area  with  weak  or 
negligible  baroclinity.  The  second  type  is  the  Pacific 
polar  low.  This  disturbance  exhibits  wavelengths  of  1000- 
1500  km  and  is  often  a  deep  system,  occasionally  exhibiting 
height  perturbations  at  300  mb  larger  than  those  at  the 
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ev  =  virtual  potential  temperature 

eQ  =  water  temperature 

uj  =  vertical  p-velocity  =  dp/dt 

Subscripts 

B  =  bottom  of  layer 
c  =  cloud  value;  center  of  grid 
m  =  mean  PBL  value 

n  =  value  for  neutral  PBL 

r  =  reference  state 

s  =  surface  value 

T  =  top  of  layer 

o  =  central  latitude;  45 N 


bulk  Richardson  number 
critical  bulk  Richardson  number 
t  i  me 

temperature 

standard  atmosphere  temperature 

component  of  wind  in  x-direction 

friction  velocity 

component  of  wind  in  y-direction 

two  dimensional  velocity  vector  =  u  i  +  v  j 

entrainment  velocity 

east-west  coordinate 

east-west  domain  extent 

north-south  coordinate 

north-south  domain  extent 

roughness  length 

inflow  angle 

df/dy 

two  dimensional  del  operator  =  a/ ax  i  +  a/ay  j 

geopotent i al 

R/Cp 

density 

terrain  pressure 
vertical  coordinate  =  p/w 
vertical  a-velocity  =  da/dt 
potential  temperature 


LIST  OF  SYMBOLS 


CQ  =  drag  coefficient 

Cy  =  friction  coefficient 

Cun  =  neutral  value  of  friction  coefficient 

C0  =  heat  transport  coefficient 

C0n  =  neutral  value  of  heat  transport  coefficient 
d-|  =  depth  of  layer  of  negative  heat  flux  at  top  of  PBL 

f  =  Coriolis  parameter 

Fu  =  friction  force  in  x-direction 

Fv  =  friction  force  in  y-direction 

g  =  gravitational  acceleration  =  9.8ms 

h  =  height  of  PBL 

he  =  equilibrium  height  of  PBL 

k  =  von  K arm an  constant  =  0.35 

Kh  =  horizontal  eddy  diffusion  coefficient 

lx  =  east-west  domain  extent 

ly  =  north-south  domain  extent 

=  moisture  convergence 
p  =  pressure 

q  =  mixing  ratio 

q*  =  saturation  mixing  ratio 

qstd  =  standard  atmosphere  saturation  mixing  ratio 

Q*  =  diabatic  heating  rate 

R  =  specific  gas  constant  for  air  =  287.04  J  kg-1  K-l 
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Time-height  cross-section  of  GE  for  G3D.  Contour 
interval  is  10  x  1 0 -  3  w  m-2.  Format  of  contours 
as  in  Fig.  C-ll  - 

Time-height  cross-sections  of  HCA  for  (a)  G8A 
and  (b)  G8D.  Contour  interval  is  20  x  1 0  "  3 
W  m“2.  Format  of  contours  as  in  Fig.  C-ll  - 

Time-height  cross-secti  ns  of  VCA  for  (aj3G8A 
and  ^b)  G8D.  Contour  interval  is  20  x  10 
W  m-2.  Format  of  contours  as  in  Fig.  C-ll  - 


C-5.  East-west  cross-sections  of  co  (dashed  line; 

heavy  solid  line  is  zero  contour),  T  (light 
solid  line)  and  trough  position  (heavy  dashed 
line)  for  (a)  G8A  and  (b)  G8D  at  84  h.  Contour 
interval  forcuis  2  ybar/s  in  (a)  and  5  ybar/s 
in  (b).  Contour  interval  for  T  is  10  K.  200  km 
of  right  side  of  domain  is  the  same  as  leftmost 


200  km  due  to  periodicity  -  155 

C-6.  12  h  accumulation  of  (a)  total,  (b)  convective 

and  (c)  stable  precipitation  in  a  fixed-frame 
reference  ending  at  108  h  for  G8D.  Contour 
interval  is  5  mm.  Grid  placement  as  in  Fig.  C-2-  158 


C-7.  Fields  of  height  and  temperature  at  (a)  850  mb, 
(b)  700  mb  and  (c)  500  mb  at  108  h  for  G8A. 
Contour  interval  for  height  is  60  in  in  (a)  and 
(b)  and  120  m  in  (c).  Grid  placement  as  in 


Fig.  C-2  . .  159 

C-8.  Similar  to  Fig.  C-7  except  for  G8D  -  160 

C-9.  Fields  of  (a)  sensible  heat  flux  and  (b)  latent 
heat  flux  at  108  h  for  G8D.  Contour  interval  is 
200  W  m  “  2 .  Grid  placement  as  in  Fig.  C-2  -  162 

C-10.  Time-height  cross-sections  of  KE  for  (a)  G8A  and 
(b)  G8D.  Contour  interval  is  10  x  1 0  3  J  m-2  in 
(a)  and  20  x  103  J  m"2  in  (b)  -  165 


C - 1 1 .  Time-height  cross-sections  of  CE  for  (a)  G8A  and 
(b)  G8D.  Contour  interval  is  20  x  10-2  w  m~2. 

Heavy  solid  line  is  zero  contour,  light  solid 

line  is  positive  contour  and  dashed  line  is 

negative  contour  -  167 

C - 1 2 .  Time-height  cross  sections  of  HCK  for  (a)  G8A  ~ 
and  (b)  G80.  Contour  interval  is  10  x  1 0 ~ 2  w  m 
in  (a)  and  20  x  1 0 ~ 2  w  m-2  in  (b).  Format  of 


contours  as  in  Fig  C-ll  -  169 

C - 1 3 .  Time-height  cross-sections  of  VCK  for  (a)  G8A 
and  (b)  G8D.  Contour  interval  is  10  x  10-2  w  m-2. 
Format  of  contours  as  in  Fig.  C-ll  -  171 

C-14.  Time-height  cross-sections  of  AE  for  (a)  G8A  and 
(b)  G8D.  Contour  interval  is  5  x  10-3  j  m-2  in 
(a)  and  10  x  10-3  j  m-2  in  (b)  -  174 
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eliminate  any  small  perturbations  which  may  have  been 
introduced  in  the  tabulation  or  interpolation  of  the  data. 
Cross-sections  of  the  mean  winter  and  genesis  jets  are  shown 
in  Fig.  2 -  1 . 

2.  Initial i zation 

The  equations  used  in  the  initialization  are  greatly 
simplified  since  initially,  all  x-deri vati ves  are  zero  and 
v=0  over  the  entire  domain.  The  initial  north-south  surface 
pressure  gradient  is  computed  geostrophi cal ly  from  the  v- 
equation  of  motion  in  a  coordinates, 

3 1  n  7T  f  u  ,  o  o  N 

3y  ~  "  RT  [d‘dl 

subject  to  the  boundary  condition  that  the  gradient  of  n  is 
zero  across  the  boundaries.  A  constant  surface  temperature 
of  282  K  is  used  and  the  1000  mb  u-components  in  Fig.  2-1 
are  specified.  The  areal  mean  surface  pressure  is  specified 
as  1013.25  mb.  Initial  geopotentials  on  pressure  surfaces 
are  computed  geostrophical ly  using 

i2-3> 

subject  to  the  boundary  condition  that  the  gradient  of  6  is 
zero  across  the  boundaries.  The  mean  height  of  each  level 
is  computed  as  follows.  Temperatures  for  both  the  mean 
winter  and  genesis  cases,  as  given  by  Mullen  (1979),  are 
tabulated  and  interpolated  to  grid  points  in  the  manner 
described  in  the  previous  section  for  the  zonal  wind. 
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wind  (m/s)  for  (a)  mean  winter  and  (b)  genesis 
jets.  Contour  interval  is  5  m/s. 


Heights  are  computed  hydrostatically  and  the  areal  mean  is 
found  by  averaging  over  each  level. 

The  initial  fields  are  then  interpolated  linearly 
with  respect  to  p0.2051  to  the  model  o  coordinates.  When 
the  model  is  integrated  with  these  fields,  a  north-south 
pressure  oscillation,  or  "sloshing"  of  3-4  mb  is  observed. 
To  achieve  a  better  balance  in  the  initial  fields,  new 
height  fields  are  computed' using  the  v-equation  of  motion 
in  a  coordinates, 


9  <t> 

3y 


fu-RT 


9  1  n  it 

ay 


(2.4) 


subject  to  the  boundary  condition  that  the  gradient  of  $  is 
zero  across  the  boundaries.  Eq.  (2.4)  must  be  solved  by  an 
iterative  process,  since  the  temperature  appears  on  the 
right  side.  On  the  first  pass,  the  temperatures  obtained 
from  Eq.  (2.3)  are  used.  On  each  succeeding  pass,  tempera¬ 
tures  obtained  from  the  previous  iteration  are  used.  The 
procedure  is  repeated  for  five  iterations.  By  balancing 
the  fields  directly  on  the  model  a  coordinates,  the  pressure 
sloshing  was  reduced  to  less  than  0.1  mb.  After  the 
initialization,  the  genesis  fields  exhibit  superad i abat i c 
lapse  rates  since  no  formal  constraint  is  imposed.  For  the 
genesis  case,  the  model  is  integrated  for  12  h  neglecting 
friction  and  diabatic  effects  to  achieve  a  balanced  state 
without  superadiabatic  layers  present.  Cross-sections  of 
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the  initialized  temperatures  for  the  mean  winter  and  genesis 
cases  are  shown  in  Fig.  2-2. 

The  final  prognostic  variable  to  be  initialized  is 
the  moisture.  For  both  jet  profiles,  the  zonal  mean  winter 
relative  humidity  fields  for  the  northern  hemisphere  from 
London  (1957)  are  used.  Values  are  tabulated  and  interpo¬ 
lated  in  the  same  manner  as  for  the  zonal  wind.  The  1000  mb 
relative  humidity  given  by  Mullen  (1979)  for  the  genesis 
stage  of  the  polar  lows  is  used.  A  cross-section  of  the 
initial  relative  humidity  is  shown  in  Fig.  2-3. 

The  remaining  variable  to  be  specified  is  the 
sea-surface  temperature  (SST).  All  experiments  apply  for 
over-water  conditions.  The  specific  SST  distribution  is 
described  with  each  experiment.  There  is  no  east-west 
gradient  in  the  SST  and  the  values  are  fixed  throughout  the 
forecast.  The  surface  moisture  is  the  saturation  value 
corresponding  to  the  SST. 

B.  ENERGY  BUDGET 

An  energy  budget  is  useful  for  studying  the  development 
and  maintenance  of  the  atmospheric  disturbance.  It  is 
common  practice  to  decompose  the  energy  equations  into 
equations  for  zonal  and  eddy  components  of  the  flow.  This 
is  particularly  useful  in  this  study,  as  a  perturbation 
grows  in  a  zonal  flow.  Energy  transfers  within  the  eddy, 
and  interactions  between  the  eddy  and  the  zonal  flow,  will 


temperature  (K)  for  (a)  mean  winter  and  (b) 
genesis  jets.  Contour  interval  is  5  K. 


PRESSURE  (MB) 


Fig.  2-3.  Initial  north-south  cross-section  of  relative 
humidity  (%).  Contour  interval  is  10%. 
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be  studied.  The  effects  of  friction  and  diabatic  processes 
also  will  be  described  in  terms  of  the  energy  budgets. 

Many  earlier  studies  have  used  the  energy  equations 
developed  by  Lorenz  (1955),  which  include  an  approximate 
form  of  the  available  potential  energy.  Gall  and  Johnson 
(1971)  demonstrated  that  this  form  of  the  equations  can 
produce  erroneous  generation  of  eddy  available  potential 
energy.  For  this  reason,  an  exact  formulation  of  the  energy 
equations,  as  derived  by  Robertson  and  Smith  (1983),  is  used 
in  the  study.  A  description  of  the  budget  equations  is 
given  in  Appendix  B. 

The  energy  budget  terms  are  calculated  using  output 
fields  from  the  forecast  model.  All  fields  are  saved  each 
6  h  and  are  t i me- a v eraged  during  the  model  integration  for 
a  one-hour  period  centered  on  the  output  time.  Therefore, 
a  12  h  forecast  field  represents  the  average  of  all  forecast 
times  between  11.5  and  12.5  hours.  All  fields  are  inter¬ 
polated  linearly  with  respect  to  In  p  from  the  model  a 
coordinates  to  constant  pressure  surfaces.  There  are  22 
constant  pressure  surfaces  in  50  mb  increments  starting  at 
25  mb  and  extending  down  to  10  75  mb.  Only  those  points  at 
or  above  the  surface  are  used  in  the  budget  calculations. 
That  is,  any  mean  quantity  at  1025  mb  is  calculated  using 
only  those  points  for  which  the  surface  pressure  >_  1025  mb. 
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The  transformation  of  the  vertical  velocity  from  the 
model  a  coordinates  to  pressure  surfaces  necessitated  the 
saving  of  several  fields.  This  t ra n s f orm at i on  is  obtained 
by  substituting  Eq.  (A-l)  into  o>  =  dp/dt,  which  yields, 

cj  =  ira  +  a  +  •  Vir  (2.5) 

O  t 

The  first  term  on  the  right  side  of  Eq.  (2.5)  is  simply 
the  model  vertical  motion  at  the  interfaces  of  the  model 
layers.  These  fields  are  interpolated  directly  to  the 
pressure  surfaces.  The  last  two  terms  are  calculated  at  the 
middle  of  each  model  a  layer  and  are  interpolated  directly 
from  these  levels  to  the  pressure  surfaces.  As  a  check,  w 
was  also  calculated  kinematically  using  the  ti me- averaged 
velocity  fields  on  pressure  surfaces.  Differences  in  the 
budget  results  were  less  than  a  few  percent. 

The  computation  of  the  reference  pressure  is  done  as 
follows.  First,  In  p  is  interpolated  linearly  with  respect 
to  9  from  the  model  a  coordinates  to  isentropic  surfaces 
at  1  K  intervals.  All  isentropic  surfaces  below  the  earth's 
surface  are  assigned  a  value  of  In  ti.  Second,  the  areal 
mean  pressure  (i.e.,  the  reference  pressure,  p  )  is 
calculated  on  each  isentropic  surface.  Third,  In  pr  is 
interpolated  to  pressure  surfaces  linearly  with  respect  to 
9.  Cross-sections  of  the  initial  efficiency  factor  for  both 
jet  structures  are  shown  in  Fig.  2-4. 
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No  horizontal  interpolation  is  necessary  for  the  model 
output  fields.  All  horizontal  derivatives  are  computed  in  a 
manner  consistent  with  the  model  finite  differencing  on  the 
scheme  C  grid.  Centered  differencing  is  used  for  vertical 
derivatives.  The  time  change  in  energy  is  computed  using  a 
centered  difference  over  a  12-h  interval. 
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III.  LINEAR  GROWTH  STUDIES 


The  concept  of  using  the  perturbation  technique  to  study 
wave  growth  in  a  zonal  flow  was  introduced  by  Charney  (1947) 
and  Eady  (1949).  These  studies  indicated  that  the  most 
unstable  wavelength  for  mean  atmospheric  conditions  was 
about  4000  km.  In  general,  this  agrees  with  observations 
which  show  that  the  ex t ra t r o p i c a  1  regions  are  dominated  by 
waves  of  this  size.  However,  more  recent  studies  by  Simons 
(1972)  and  Gall  (1976a)  indicate  that  the  most  unstable 
wavelength  depends  critically  on  the  mean  flow  structure.  Gall 
(1976a)  showed  with  the  linearized  equations  that  wavenumber 
15  was  most  unstable  for  mean  winter  conditions.  However, 
wavenumber  15  failed  to  dominate  over  longer  wavelengths  due 
to  wave-zonal  flow  interaction  and  surface  friction  effects 
for  finite  ampl i tude . 

The  purpose  of  this  chapter  is  to  study  the  growth  rates 
for  the  mean  winter  and  genesis  jets  described  in  Chapter 
II.  First,  linear  growth  rates  will  be  computed  for  both 
jets  and  the  energetics  and  structure  of  the  most  unstable 
waves  will  be  examined.  Second,  the  effect  of  wave-zonal 
flow  interaction  will  be  determined.  Third,  the  effect  of 
linearized  surface  friction  on  the  growth  rates  will  be 
examined. 
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A.  EXPERIMENTAL  DESIGN 

A  reduced  version  of  the  numerical  model  described  in 
Appendix  A  is  used  to  study  the  linear  growth  rates  of 
perturbat i ons  in  the  mean  winter  and  genesis  jets.  Two 
simplifications  are  made  to  the  model  in  these  experiments. 
First,  the  physical  processes  are  not  included.  Second, 
after  each  time  step,  the  zonal  averages  of  u,  v,  T  and  u 
are  replaced  with  the  initial  zonal  averages  of  their 
respective  values,  which  eliminates  interactions  between  the 
growing  waves  and  the  zonal  flow.  These  simplifications 
transform  the  nonlinear  model  into  a  linear  model  with 
respect  to  interaction  with  the  mean  flow.  The  nonlinear 
interactions  between  waves  are  still  included. 

The  Cartesian  grid  has  a  horizontal  resolution  of  100  km 
with  94  grid  intervals  in  the  x-direction  and  64  intervals 
in  the  y-direction.  Therefore,  the  east-west  extent  of  the 
grid  is  9400  km  or  1/3  of  the  global  circumference  at  46N. 
Since  the  genesis  jet  exhibited  the  largest  growth  rates  for 
the  shortest  waves,  another  version  of  the  model  was  run  for 
the  genesis  case  only.  This  version  uses  75  km  resolution 
with  94  intervals  in  the  x-direction  and  85  intervals  in  the 
y-direction.  The  east-west  extent  of  this  grid  is  7050  km 
or  1/4  of  the  global  circumference  at  45N.  The  model 
vertical  levels  are  given  in  Table  A-l. 
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Initial  v-component  perturbat i ons  are  introduced  for 
model  wavenumbers  1-10.  These  correspond  to  earth  wave- 
numbers  3-30  (4-40)  in  increments  of  3(4)  for  the  100(75)  km 
grid.  The  perturbat i on ,  which  is  constant  in  the  vertical, 
is  given  by , 


=  V  sin  (^p)sin£[ 


(y-y„  )* 


-] 


(3.1) 


The  y-structure  is  chosen  so  v  '  =  0  on  the  north  and  south 
boundaries.  The  initial  amplitude  of  the  perturbations  is 
0.1  m  s-1  for  the  mean  winter  jet.  However,  the  initial 
amplitude  for  the  genesis  jet  was  set  to  0.01  m  s-1  because 
it  exhibited  such  high  growth  rates.  This  insured  a  longer 
period  over  which  to  compute  the  linear  growth  rate.  As 
noted  by  Gall  (1976a),  the  imposition  of  the  initial  pertur¬ 
bation  leads  to  generation  of  gravity  waves.  Due  to  the 
presence  of  these  waves,  it  was  estimated  that  the  growth 
rate  was  computed  to  within  an  accuracy  of  0.05  day-1  in 
this  study . 

The  technique  used  to  evaluate  the  linear  growth  rate  is 
similar  to  that  used  by  Gall  (1976b),  who  inserted  only  one 
wave  per  forecast.  Therefore,  Gall  required  many  forecasts 
to  achieve  the  growth  rate  spectrum  fo*-  ttie  given  zonal 
flow.  In  this  experiment,  10  wavenumbers  are  included 
initially,  so  that  interactions  between  waves  may  alter  the 
growth  rates  of  some  of  the  waves.  However,  these  nonlinear 
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interactions  will  remain  small  since  the  pe r t u r ba t i o n s  are 
small.  As  a  check,  the  model  was  run  with  only  one  wave- 
number  in  the  initial  fields  for  two  different  wavenumbers. 
The  growth  rates  computed  from  these  runs  were  identical  to 
those  computed  when  all  10  wavenumbers  were  included,  which 
verified  that  the  nonlinear  effects  were  small. 

The  growth  rates  and  ph<-se  speeds  were  computed  as 
follows.  A  double  Fourier  decomposition  of  the  v-component 
in  the  lowest  model  layer  was  calculated  every  12  h  to  get 
the  amplitude  and  phase  of  each  of  the  10  waves.  During  the 
period  of  linear  growth,  the  growth  rate  for  each  wave  was 
computed  as 


v 


3  ( 1  n  v ') 
3  t 


Inv1  -  lnv' 
e _ s 

At 


(3.2) 


where  the  subscripts  s  and  e  refer  to  the  start  and  the  end 
of  tne  linear  growtn  period.  Phase  speeds  were  also  calcu¬ 
lated  for  each  1 2 - h  period,  and  were  averaged  during  the 
period  of  linear  growth  fur  each  wave. 

B.  WAVE  GROWTH 

1 .  Linear  Growth  Rates  and  Wave  Structure 

The  growth  rates  computed  for  the  mean  winter  and 
genesis  jets  are  shown  in  Fig.  3-1.  There  are  two  major 
differences  between  the  curves.  First,  the  genesis  jet 
exhibits  much  larger  growth  rates,  especially  for 
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previous  experiments,  the  physical  processes  are  neglected 
and  the  zonal  averages  of  u,  v,  T,  and  u  are  restored  eacn 
time  step . 

The  time  evolution  of  the  terms  CE,  HCK  and  VCK  for 
W12,  G8  and  G40  are  given  in  Figs.  3-7,  3-3  and  3-9,  respec¬ 
tively.  For  W12,  CE  is  the  dominant  term,  which  shows  that 
the  growth  of  this  wave  is  primarily  through  conventional 
baroclinic  instability.  HCK  and  VCK  are  5-6  times  smaller 
than  CE,  which  indicates  that  barotropic  and  vertical  shear 
effects  are  small.  Also,  HCK  and  VCK  tend  to  cancel  since 
they  are  nearly  the  same  size  but  of  opposite  sign.  For  G3 , 
CE  is  again  the  dominant  term.  VCK  is  positive  and  about 
one-third  of  CE  ,  which  suggests  that  Kel v i n-Hel mhol tz 
instability  may  play  a  secondary  role  is  the  development  of 
this  disturbance.  The  HCK  is  again  negative,  but  is  only 
about  half  the  magnitude  of  VCK.  For  G40,  both  CE  and  VCK 
are  positive  and  nearly  equal  in  magnitude,  which  shows  that 
this  wave  draws  as  much  energy  from  KZ  as  from  AE.  Term  HCK 
is  negative  and  an  order  of  magnitude  less  than  CE  and  VCK. 

It  is  apparent  from  Figs.  3-7,  3-8  and  3-9  that  tne 
length  of  time  of  linear  growth  of  G40  is  significantly  less 
than  that  of  the  longer  waves  G8  and  W12.  While  the  linear 
growth  rates  extended  past  120  h  for  G8  and  W12,  they  only 
continued  to  about  60  h  for  G40.  The  reasons  for  the  rapid 
stabilization  of  this  wave  will  be  discussed  in  the 
fol 1 owi ng  section. 


3 u / ay  <  0,  a  positive  (negative)  correlation  between  u*  and 
v'  leads  to  an  increase  (decrease)  in  KE.  Generally, 
troughs  which  tilt  opposite  the  horizontal  wind  shear  are 
are  barotropi cal ly  unstable.  VCK  represents  the  term 
( -u 1  a)  ‘  3u/ ap)  and  therefore  represents  the  same  type  of 
motion  as  HCK  but  in  the  vertical  plane.  If  the  wind 
increases  with  height,  a  positive  (negative)  correlation 
between  u'  and  u*  will  lead  to  an  increase  (decrease)  in  KE 
This  is  referred  to  as  Kel v i n-He 1 mhol t z  (shearing)  instabil 
ity  which  can  become  significant  for  R.  <  1  (Orlanski, 
1968).  The  term  VCK  is  not  included  in  q ua s i - g eos t r oph i c 
models  such  as  the  one  used  by  Sardie  and  Warner  (1983). 

The  most  unstable  mode  for  the  mean  winter  jet  is 
wavenumber  12,  referred  to  as  W12,  which  corresponds  to  a 
wavelength  of  about  2400  km.  For  the  genesis  jet,  the  most 
unstable  waves  are  wavenumbers  8  (G8)  and  40  (G40),  which 
correspond  to  wavelengths  of  about  3500  and  720  km, 
respectively.  Three  different  model  domains  are  used  to 
separately  simulate  these  waves  (Table  3-1).  As  in  the 

Table  3-1.  Number  of  grid  intervals  and  horizontal 
resolution  for  waves  W12,  G8  and  G40. 


Wave 

x- i nterval s 

y-  i  nterval s 

Hori zontal 
Resol uti on  (km) 

W12 

24 

64 

100 

G8 

35 

64 

100 

G40 

18 

160 

40 

50 


speeds  in  the  genesis  case.  For  the  mean  winter  jet,  the 
maximum  phase  speed  occurs  at  wavenumber  9,  and  there  is  a 
general  decrease  in  phase  speed  for  shorter  waves.  For  the 
genesis  jet,  the  maximum  phase  speed  is  also  at  wavenumber  9 
with  a  similar  decrease  in  phase  speed  for  the  shorter 
waves.  The  lower  phase  speeds  for  the  shorter  waves  are 
because  these  shallow  systems  have  a  lower  steering  level. 

2 .  Energy  Conversions 

The  most  unstable  mode  for  the  mean  winter  jet  and 
the  two  most  unstable  modes  for  the  genesis  jet  are  now 
examined  in  more  detail.  For  each  of  these  waves,  the  model 
is  configured  to  allow  just  that  wave  in  the  domain,  so  the 
process(es)  leading  to  its  development  can  be  isolated.  In 
particular,  the  baroclinic  conversion  (CE)  between  eddy 
available  potential  energy  (AE)  and  eddy  kinetic  energy  (KE) 
and  the  conversion  (CK)  between  zonal  kinetic  energy  (KZ) 
and  KE  are  examined. 

The  baroclinic  term,  CE,  is  the  correlation  between 
vertical  velocity  and  temperature  perturbations,  ( u> 1  T  ‘  )  .  KE 
increases  for  rising  warm  air  and  sinking  cold  air  at  the 
same  elevation,  and  vice  versa.  The  conversion  CK  is 
separated  into  the  horizontal  (HCK)  and  vertical  (VCK) 
components.  Since  v=0  in  this  linear  study,  HCK  represents 
the  barotropic  term  (u'v'3u/3y)  which  is  the  correlation  of 
u  and  v  component  perturbations  multiplied  by  the  meridional 
gradient  of  the  mean  zonal  wind.  North  of  the  jet,  where 
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4.  Phase  angles  for  wavenumbers  (a)  6,  (b)  12, 

(c)  18,  (d)  24,  and  (e)  30  for  the  mean  winter 


wavelength  of  maximum  growth  to  shorter  wavelengths.  The 
differences  between  the  growth  rates  presented  here  for 
the  mean  winter  and  genesis  jets  is  evidence  of  this. 

Staley  and  Gall  (1977)  also  found  that  wavelengths  less 
than  3000  km  (wavenumber  >10)  are  confined  below  500  mb 
and  exhibit  their  maximum  amplitude  at  the  surface. 

The  phase  structure  of  wavenumbers  6,  12,  18,  24  and 

30  for  the  mean  winter  jet  and  wavenumbers  6,  12,  18,  24,  30 

and  40  for  the  genesis  jet  are  shown  in  Figs.  3-4  and  3-5. 

For  both  jets,  wavenumber  6  has  a  pronounced  tilt  to  the 
west  with  height  between  the  surface  and  700  mb  and  a 
reduced  westward  tilt  up  to  300  mb.  For  wavenumbers  12  and 
18  in  the  mean  winter  case  and  wavenumber  12  in  the  genesis 
case,  a  moderate  westward  tilt  is  observed  in  the  lowest 
400  mb  with  little  or  no  tilt  -above.  For  the  remaining 
shorter  waves,  a  very  strong  westward  tilt  is  found  in  the 
lowest  200-400  mb.  This  tilt  becomes  more  pronounced  for 
increasing  wavenumber.  Since  a  westward  tilt  is  associated 
with  baroclinic  processes,  these  results  suggest  that  baro- 
clinic  instability  is  exhibited  by  all  waves,  but  is 
confined  to  the  lowest  few  hundred  millibars  for  wavenumbers 
greater  than  12. 

The  phase  speeds  for  the  perturbations  in  the  mean 
winter  and  genesis  jets  are  shown  in  Fig.  3-6.  Phase  speeds 
in  the  genesis  jet  case  are  30-40%  higher  than  those  in  the 
mean  winter  jet  as  would  be  expected  from  the  larger  wind 


at  60  h  for  (bj  genesis  jet.  Amplitudes  nave  ueen 
normalized  so  that  the  largest  amplitude  for  a  given 
wavenumber  is  100.  Contour  interval  is  20. 


that  small,  low-level  static  stability  and/or  an  increase  in 
the  low-level  wind  shear  (i.e.,  a  small  R j )  tends  to  shift 
the  maximum  growth  rate  to  shorter  wavelengths.  It  can  be 
seen  from  Fig.  3-2  that  the  mean  winter  case  has  <  10 
over  a  large  depth  of  the  troposphere  in  the  vicinity  of  the 
jet  and  in  the  low  levels  poleward  of  the  jet.  The  genesis 
jet  case  exhibits  the  same  pattern  with  Rn-  <  2  in  the  low- 
levels  poleward  of  the  jet  core.  In  fact,  many  values  are 
less  than  1.  This  difference  in  R .  is  consistent  with  the 
large  difference  in  growth  rates  between  the  mean  winter  and 
genesis  jet  cases  for  wavenumbers  greater  than  21. 

The  normalized  amplitude  of  the  waves  during  the 
period  of  linear  growth  is  similar  for  the  mean  winter  and 
the  genesis  jet  cases  (Fig.  3-3).  Wavenumber  6  has  a  maxi¬ 
mum  amplitude  in  the  lowest  layer  and  also  near  300  mb.  The 
shorter  waves  (wavenumber  >  15)  are  confined  to  the  lower 
troposphere  with  the  maximum  amplitude  in  the  lowest  layer. 
The  amplitude  of  these  short  waves  decreases  to  20%  of  the 
near-surface  values  at  57  0  mb  for  the  mean  winter  jet  and  at 
650  mb  for  the  genesis  jet. 

These  results  are  consistent  with  those  of  Staley 
and  Gall  (  1977  ),  who  used  a  4-level  qua s i - geos t ro ph i c  model 
to  study  the  effect  of  wind  shear  and  static  stability  on 
baroclinic  development.  They  found  that  decreasing  the  low- 
level  static  stability  (i.e.,  decreasing  R-j  )  shifts  the 


wavenumbers  greater  than  20.  Second,  the  most  unstable 
wavelength  is  markedly  different  for  the  two  jets.  For  the 
mean  winter  jet,  wavenumber  12  is  most  unstable,  with  a 
growth  rate  of  0.59  day-1.  For  the  genesis  jet,  wavenumber 
40  is  most  unstable  with  a  growth  rate  of  1.32  day-1,  and  a 
secondary  maximum  of  1.29  day-1  at  wavenumber  8.  Wave- 
numbers  larger  than  40,  which  were  not  tested  here,  may 
exhibit  even  larger  growth  rates. 

The  growth  rates  for  the  genesis  jet  were  relatively 
easy  to  determine.  However,  growth  rates  for  the  mean 
winter  jet  for  wavenumbers  greater  than  24  were  more  diffi¬ 
cult  to  estimate.  These  waves  exhibited  little  or  no  growth 
for  up  to  7  days  and  then  experienced  rapid  growth.  The 
longer  waves  (wavenumbers  £  24)  generally  began  exponential 
growth  between  days  2  and  3  of  the  forecasts.  It  was  judged 
that  the  eventual  growth  of  the  shortest  waves  was  due  to 
nonlinear  interactions  with  the  longer  waves,  which  by  day  7 
had  reached  amplitudes  >  5  m  $-1.  Model  runs  including  only 
wavenumber  28,  or  including  only  wavenumber  40  verified  that 
the  linear  growth  of  these  waves  are  quite  small  for  the 
mean  winter  jet. 

One  explanation  for  the  shorter  wavelength  disturb¬ 
ance  in  the  genesis  jet  is  illustrated  by  the  Richardson 
number  (R..)  distributions  for  both  jets  (Fig.  3-2).  R..  is 

the  ratio  of  the  static  stability  to  the  square  of  the 
vertical  wind  shear.  Staley  and  Gall  (1977)  demonstrated 


or  set  A  experiments. 


ENERGY  CONVERSION  IWM“ 


East-west  cross  sections  of  vertical  velocity, 
temperature  and  trough  position  for  W12,  G8  and  G40  are 
given  in  Figs.  3-10,  3-11  and  3-12,  respectively.  Very 


similar  structures  are  observed  for  W12  and  G8.  For  both  of 
these  waves,  the  trough  tilts  to  the  west  with  height  with 
the  maximum  slope  from  800  to  600  mb.  The  warmest  air  is 
found  ahead  of  the  trough  with  a  slight  tilt  to  the  east 
with  height  below  700  mb  and  a  slight  tilt  to  the  west  with 
height  above  700  mb.  The  maximum  upward  motion  is  found  at 
750  mb.  A  broad  area  of  relatively  weak  subsidence  occurs 
behind  the  W12  and  G8  troughs. 

The  structure  of  the  G40  wave  is  significantly 
different  from  W 1 2  and  G8  in  that  the  warmest  air  slopes  to 
the  west  with  height,  as  does  the  trough.  Therefore,  only  a 
shallow  layer  (about  250  mb  thick)  of  warm  air  is  found 
above  the  position  of  the  surface  low  with  relatively  cold 
air  above  750  mb.  The  upward  motion  above  the  position  of 
the  surface  low  leads  to  the  generation  of  KE  below  750  mb 
(i.e.,  CE  >  0)  and  a  destruction  of  KE  above  750  mb.  This  is 
consistent  with  the  shallow  nature  of  this  wave. 

The  relatively  large  values  of  VCK  for  G8  and  espe¬ 
cially  for  G40  are  due  to  the  strong  vertical  wind  snear  in 
the  genesis  jet.  Since  VCK>0  and  8u/8p<0,  this  implies  a 
positive  correlation  exists  between  u1  and  w 1  .  Fields  of  u' 


55 


Fig.  3-10.  East-west  cross-section  of  w  (dashed  line;  heavy 
dark  line  is  zero  contour),  T  (solid  line)  and 
trough  position  (heavy  dashed  line)  for  wave- 
number  12  in  mean  winter  jet.  Contour  interval 
is  1  ybar/s  for  u  and  10  K  for  T.  200  km  of 
right  side  of  domain  is  the  same  as  leftmost 
200  km  due  to  periodicity. 
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(KM/100) 


0  except  for  wavenumber  8 
itour  interval  is  2  ubar/s 


PRESSURE  (MB) 


Fig.  3-12.  Similar  to  Fig.  3-10  except  for  wavenumber  40 
in  genesis  jet.  Contour  interval  is  5  ubar/s 
for  oj  and  10  K  for  T.  80  km  of  right  side  of 
domain  is  the  same  as  leftmost  80  km  due  to 
periodicity. 


58 


and  w 1  for  G40  (Fig.  3-13)  exhibit  a  pattern  of  upward 
motion  in  areas  of  east  winds  and  downward  motion  in  areas 
of  west  winds. 

Orlanski  (1968)  showed  that  Kel v i n-Hel mhol tz  insta¬ 
bility  can  dominate  over  conventional  baroclinic  instability 
for  small  .  Stone  (1966)  showed  that  helical  circulations 
parallel  to  the  zonal  flow  dominate  for  0.95  <  R^<  0.25 
(symmetric  instability),  whereas  the  circulations  are 
perpendicular  to  the  zonal  flow  for  R-j  <  0.25  (Kelvin- 
Helmholtz  instability).  A  north-south  cross-section  of  v 
and  w through  G40  is  given  in  Fig.  3-14.  Strong  vertical 
shear  in  the  v-component  coupled  with  the  horizontal 
gradient  of  u  leads  to  a  circulation  in  the  vertical  plane 
about  a  horizontal  axis  near  800  mb.  Since  this  is  a  north- 
south  cross-section,  this  circulation  is  parallel  to  the 
zonal  flow  and  suggests  that  symmetric  instability  is 
present.  However,  the  u1  and  w'  fields  in  G40  (Fig.  3-13) 
suggest  that  this  disturbance  is  strongly  tilted  by  the 
horizontal  wind  shear  and  therefore  this  helical  circulation 
may  also  be  driven  by  the  Kel v i n-Hel mhol tz  instability.  R^ 
values  computed  in  the  lowest  level  range  from  0.03-2.0, 
which  indicates  that  either  type  of  motion  is  possible.  At 
any  rate,  the  presence  of  this  helical  circulation  is 
consistent  with  the  results  of  Stone  (1966)  and  Orlanski 
(1968).  Therefore,  based  on  the  energy  conversions  and  wave 
structure,  it  is  concluded  that  G40  is  a  "mixed"  baroclinic 
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subsection  of  the  forecast  domain.  Southern  boundary  of  figures  is 
at  the  center  of  the  forecast  domain.  80  km  of  right  side  of  domain 
is  the  same  as  the  leftmost  80  km  due  to  periodicity. 
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wave.  That  is,  energy  is  gained  from  both  KZ  and  AE. 
Although  VCK  is  part  of  the  barotropic  term,  it  is 
considered  a  baroclinic  process  since  it  involves  the 
vertical  wind  shear. 

C.  EFFECT  OF  WAVE-ZONAL  FLOW  INTERACTION 

Although  Gall  (1976a)  showed  that  the  most  unstable 
growth  rate  occurred  at  wavenumbers  12-15,  the  dominant 
scale  of  motion  was  eventually  found  to  be  wavenumbers  5-7. 
Subsequestly ,  Gall  (  1976b)  included  wave-zonal  flow  inter¬ 
action  and  surface  friction  effects.  It  was  found  that  the 
growth  rate  in  the  lower  troposphere  exhibited  the  most 
rapid  reduction  due  to  increasing  static  stability  and 
reduced  vertical  shear.  Coupled  with  the  effects  of  surface 
friction  in  the  lower  troposphere,  short  waves  (wavenumbers 
12-15)  stopped  growing  earlier  than  longer  waves  which 
extended  throughout  the  troposphere.  Gall  (1976b)  concluded 
that  wave  interaction  with  the  zonal  flow  plus  surface 
friction  modified  the  waves  and  the  zonal  flow  so  that 
wavenumbers  5-7  dominated  the  flow. 

It  was  observed  in  the  experiments  with  the  genesis  jet 
in  the  previous  section,  that  after  G40  had  developed,  other 
disturbances  somewhat  larger  than  wavenumber  40  developed. 

It  was  suspected  that  a  maximum  in  the  growth  rate  spectrum 
existed  between  wavenumber  8  and  40  for  a  slightly  modified 
genesis  jet.  Therefore,  in  this  section,  the  effects  of 


wave-zonal  flow  interaction  are  studied  for  the  mean  winter 
and  genesis  jets.  These  experiments  (set  B)  differ  from  the 
previous  experiments  (set  A)  in  the  following  way.  First, 
the  model  is  integrated  in  time  until  the  linear  growth  rate 
of  the  waves  ceases.  This  occurs  at  120  h  for  the  mean 
winter  jet  and  at  84  h  for  the  genesis  jet.  Second,  the 
zonal  means  are  computed  for  u,  v,  T  and  ir  from  these  fore¬ 
cast  fields  and  are  used  as  the  initial  fields.  The  growth 
of  the  waves  in  the  modified  mean  winter  and  genesis  jets  is 
then  calculated  exactly  as  in  the  set  A  experiments. 

The  growth  rates  for  the  set  B  experiments  are  compared 
in  Fig.  3-15  with  those  from  set  A.  For  the  mean  winter  jet 
modified  by  the  wave-zonal  flow  interaction  the  growth  rates 
are  reduced  by  10-30%  for  wavenumbers  greater  than  9.  There 
is  a  tendency  for  the  growth  rates  of  the  highest  wave- 
numbers  to  have  the  greatest  reduction.  Very  significant 
changes  occurred  with  the  modified  genesis  jet.  Here  again, 
growth  rates  are  reduced  for  all  wavenumbers  greater  than  9 
with  the  most  significant  reduction  at  the  shortest  wave¬ 
lengths.  In  fact,  the  growth  rate  for  wavenumber  40 
decreased  from  1.82  day“l  to  0.82  day-1.  In  the  previous 
section  it  was  shown  that  this  wave  involves  both  baroclinic 
instability  and  Kel v i n-He 1 mhol tz  instability.  Stone  (1966) 
and  Orlanski  (1968)  demonstrated  that  must  be  <  1  for  the 
latter  to  exist  and  this  condition  was  satisfied  (cf,  Fig. 
3-2)  for  the  set  A  experiments.  However,  the  modified 
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genesis  jet  for  the  set  B  experiment  does  not  meet  this 
condition.  R  values  in  the  v i c i n ty  of  the  modified 
genesis  jet  in  the  lowest  levels  are  in  the  range  of  2-5 
(Fig.  3-16),  which  is  a  small  but  significant  increase 
relative  to  the  set  A  experiment.  This  R .  increase  is  due 
to  an  increase  of  the  static  stability  in  the  lowest  levels. 
The  decreased  growth  rates  for  the  shortest  waves  are 
consistent  with  the  increased  R .  values  (Stone,  1966; 
Orlanski,  1968;  Staley  and  Gall,  1977).  R.  values  poleward 
of  the  jet  in  the  lowest  levels  are  unchanged  from  the  set  A 
experiments  and  remain  less  than  2. 

The  growth  rate  for  wavenumber  20  for  the  genesis  jet 
was  not  reduced  as  much  as  the  other  waves.  In  fact,  the 
growtn  rate  spectrum  for  the  modified  genesis  jet  exhibits  a 
maximum  at  wavenumber  8  (1.17  day"1)  and  also  at  wavenumber 
20  (1.2  day”1).  It  is  suggested  that  the  peak  at  wavenumber 
20  did  not  appear  in  the  A  experiment  due  to  the  magnitude 
of  the  growth  rate  for  the  "mixed"  baroclinic  mode. 

To  study  wavenumber  20  (G20),  the  model  was  initialized 
with  just  this  wave  in  the  domain.  Here,  a  horizontal 
resolution  of  50  km  with  28  intervals  in  the  x-direction 
and  128  intervals  in  the  y-direction  is  used.  The  initial 
fields  were  generated  as  in  the  B  experiments  described 
earl i er . 
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Fig.  3-16.  Initial  zona  1 1 y- a ve ra g ed  Richardson  number 

for  modified  genesis  jet.  Contour  interval  is 
2.  Values  above  10  are  not  contoured. 
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The  time  evolution  of  CE ,  HCK  and  VCK  for  G20  is  shown 
in  Fig.  3-17.  Term  CE  is  dominant,  which  indicates  that 
this  wave  receives  energy  through  conventional  baroclinic 
processes.  Term  VCK  is  also  positive  but  is  only  one-fourth 
of  CE.  Term  HCK  is  negative  and  about  one-half  the  magni¬ 
tude  of  VCK . 

An  east-west  cross-section  of  vertical  motion,  tempera¬ 
ture  and  trough  position  for  G20  is  given  in  Fig.  3-18.  The 
trough  tilts  westward  with  height  with  the  largest  slope 
between  850  and  700  mb.  The  warmest  air  extends  from  the 
surface  trough  position  to  750  mb  with  little  or  no  tilt. 
Above  750  mb,  the  warmest  air  tilts  to  the  west.  With  the 
strongest  vertical  motion  occurring  just  east  of  the  surface 
trough  position,  generation  of  KE  is  a  maximum  in  the  lowest 
250  mb.  This  pattern  is  similar  to  that  observed  in  G40  in 
the  previous  section  and  is  consistent  with  the  shallow 
nature  of  G20  and  G40  compared  to  W12  and  G8 . 

The  surface  pressure  pattern  for  G20  is  shown  in  Fig. 
3-19.  A  well-defined  surface  low  develops  about  1000  km 
north  of  the  jet  stream.  Although  the  horizontal  scale  of 
G20  is  smaller  than  that  of  ex t r a t ro p i ca 1  ,  sy nopt i c - sc  a  1  e 
waves,  surface  pressure  patterns  in  both  waves  are  similar. 
In  addition,  the  energetics  for  both  waves  are  similar  with 
large  positive  CE  and  smaller,  negative  HCK. 
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air-sea  temperature  difference  fields  are  shown  with  the 
individual  experiments.  Surface  vapor  pressures  are  the 
saturation  values  correspond i ng  to  the  sea-surface  tempera¬ 
tures.  Sea-surface  temperature  and  moisture  values  remain 
constant  over  each  forecast. 

The  diabatic  runs  use  the  complete  model  while  the 
adiabatic  runs  neglect  precipitation  physics  and  sensible 
and  latent  heating  within  the  PBL.  Also,  the  adiabatic  runs 
utilize  the  surface  momentum  fluxes  to  incorporate  surface 
friction  and  the  surface  heat  flux  to  predict  the  PBL  depth 
even  though  the  sensible  heating  due  to  the  surface  heat 
flux  is  not  added  to  the  temperature  field. 

A.  WAVE  GROWTH  AND  STRUCTURE 

A  graph  of  the  minimum  central  pressure  versus  time  of 
G20  for  both  the  adiabatic  (G20A)  and  diabatic  (G20D)  simu¬ 
lations  is  given  in  Fig.  4-1.  The  inclusion  of  diabatic 
effects  results  in  G20D  attaining  a  central  pressure  of 
984  mb,  which  is  19  mb  lower  than  G20A.  Clearly,  diabatic 
effects  serve  to  intensify  the  disturbance.  This  is 
expected  based  on  the  studies  of  Bullock  and  Johnson  (1971), 
Tracton  (1973),  Gall  (1976c)  and  Sanders  and  Gyakum  (1980). 

Mullen  (1983)  computed  deepening  rates  for  several 
Pacific  polar  lows  according  to  the  method  proposed  by 
Sanders  and  Gyakum  (1980),  who  defined  1  Bergeron  (B)  as  a 
central  pressure  fall  of  1  mb  h~l  for  24  h  at  60°  latitude. 


presented  in  Appendix  C.  The  results  of  the  experiments  in 
Appendix  C  are  used  for  comparison  of  the  polar  low  to  an 
extratropi cal  sy nopt i c - seal e  system. 

The  following  procedure  is  used  to  incorporate  the 
initial  perturbation  into  the  fields  for  the  model  simula¬ 
tions.  First,  the  perturbation  for  a  single  wave  is  imposed 
in  the  initialized  fields  in  the  manner  described  in  Chapter 
II  (i.e.,  constant  amplitude  in  the  vertical).  Second,  the 
model  is  integrated  for  36  h  with  the  initial  zonal  mean 
values  of  u,  v,  T,  q  and  n  restored  every  time  step  under 
inviscid,  adiabatic  conditions.  This  allows  the  wave  to 
reach  the  point  at  which  exponential  growth  starts  as  in  the 
Chapter  III  experiments.  This  insures  the  proper  vertical 
structure  of  the  wave  at  the  initial  time.  These  36  h 

forecast  fields  serve  as  the  initial  fields  for  all  simula¬ 

tions  presented  here. 

The  sea-surface  temperature  for  these  experiments  is 
computed  as  follows.  The  mean  air-sea  temperature  differ¬ 
ence  given  by  Mullen  (1979)  for  the  genesis  conditions  is 
interpolated  to  the  forecast  grid  in  the  same  manner  as  for 

the  u-component  described  in  Chapter  II.  Mullen  (1979) 

found  the  maximum  air-sea  temperature  difference  of  6.8°C  to 
be  approximately  750  km  north  of  the  jet  axis  and  decreased 
to  about  4.5°C  to  the  north  and  to  about  Z°C  south  of  the 
jet  axis.  This  difference  is  added  to  the  lowest  level 
initial  temperature  field.  The  sea-surface  temperature  and 
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IV.  NONLINEAR  MODEL  SIMULATIONS 


It  was  shown  in  Chapter  III  that  the  linear  growth  rates 
for  both  the  mean  winter  and  genesis  jets  are  largest  for 
relatively  large  (wavenumber  8-9)  baroclinic  waves  when 
frictional  effects  and  wave-zonal  flow  interaction  are 
considered.  In  addition,  the  genesis  jet  exhibits  large 
growth  rates  for  two  smal 1 er- seal e  systems.  The  first, 
wavenumber  40,  was  found  to  be  very  shallow  and  forms 
through  baroclinic  and  Kelvin-Helmholtz  instabilities.  After 
this  wave  modifies  the  zonal  flow  in  the  vicinity  of  the 
jet,  a  peak  in  the  growth  rate  spectrum  occurs  at  wavenumber 
20  (1400  km).  This  particular  wave,  G20,  forms  through 
conventional  baroclinic  instability  and  resembles  the 
Pacific  polar  low  described  by  Reed  (1979)  and  Mullen 
(1979).  However,  the  vertical  scale  of  G20  is  smaller  than 
observed  for  the  Pacific  polar  low. 

The  development  and  structure  of  G20  will  now  be  studied 
in  detail.  This  will  be  accomplished  through  the  examina¬ 
tion  of  model  output  and  energy  budgets  of  the  storm  at 
various  stages  of  development  for  both  adiabatic  and 
diabatic  simulations.  The  model  configuration  to  study 
these  waves  is  the  same  as  that  used  in  studying  the  energy 
transfer  in  these  waves  in  Chapter  III  (Table  3-1).  A 
similar  study  of  wavenumber  8  for  the  genesis  jet  is 
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The  second  major  point  is  that  the  linear  growth  rate 
spectrum  is  strongly  affected  by  processes  which  act  to 
continuously  modify  the  initial  zonal  flow.  The  processes 
tested  here  include  wave-zonal  flow  interation  and  surface 
friction.  It  was  found  that  surface  friction  damps  all 
wavelengths,  but  is  especially  effective  in  damping  the 
shortest  waves.  In  addition,  surface  friction  shifted  the 
most  unstable  wavenumber  for  the  mean  winter  jet  from  wave- 
number  12  to  9.  It  was  also  found  that  wave-zonal  flow 
interaction  modified  the  shortest  waves  the  most  and  led  to 
a  significant  change  in  the  growth  rate  spectrum  for  the 
genesis  jet.  The  baroclinic/Kelv in-Helmholtz  wave  was 
damped  considerably,  as  the  growth  rate  was  reduced  from 
1.82  day-1  to  0.82  day-1.  This  was  due  to  an  increase  in 
the  low-level  R.  values  from  less  than  1  to  values  of  2-5. 

i 

The  modified  genesis  jet  exhibited  maximum  growth  rates  at 
wavenumbers  8  and  20.  It  was  found  that  wavenumber  20  was  a 
baroclinic  wave  with  maximum  amplitude  in  the  lowest  level 
and  resembles  the  Pacific  polar  low  described  by  Reed  (1979) 
and  Mullen  (1979).  In  addition,  this  wave  is  not  as  deep  as 
the  longer  baroclinic  waves,  with  the  maximum  upward  motion 
occurring  at  800  mb  as  opposed  to  750  mb  for  the  longer 
waves.  Friction  is  found  to  damp  wavenumber  20  strongly, 
which  suggests  that  polar  lows  cannot  attain  observed 
intensities  through  baroclinic  instability  alone. 
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the  growth  rate  spectrum  of  the  original  genesis  jet  such 
that  the  synoptic-scale  wave  dominates.  Clearly,  the  deve¬ 
lopment  of  a  polar  low  must  involve  processes  other  than 
baroclinic  instability  to  attain  observed  intensities.  The 
effect  of  diabatic  processes  on  the  polar  low  will  be 
studied  in  Chapter  IV. 

E.  CHAPTER  SUMMARY 

The  experiments  described  in  this  chapter  lead  to  two 
important  points.  First,  significant  differences  exist  in 
the  growth  rate  spectrum  for  the  two  jet  profiles.  It  was 
found  that  the  mean  winter  jet  yields  a  maximum  growth  rate 
for  wavenumber  12  (0.59  day-1).  This  result  supports  the 
findings  of  Gall  (1976a),  who  computed  linear  growth  rates 
for  a  mean  zonal  flow  similar  to  the  mean  winter  jet.  He 
found  wavenumbers  12-15  to  be  most  unstable  with  a  growth 
rate  of  0.6  day-1.  In  addition,  the  shape  of  the  curve  in 
Fig.  3-1  for  the  mean  winter  jet  is  quite  similar  to  that 
given  by  Gall  (1976a),  who  included  only  wavenumbers  <  15. 
The  genesis  jet  case  yields  a  maximum  growth  rate  for  wave- 
number  8  (1.22  day-1)  and  40  (1.82  day-1).  It  was  shown 
that  wavenumber  12  for  the  mean  winter  case  and  wavenumber 
8  for  the  genesis  case  are  baroclinic  waves  which  extend 
throughout  the  troposphere.  However,  wavenumber  40  was 
found  to  be  a  barocl i nic/Kel vi n-Helmhol tz  wave  which  is 
confined  to  the  surface-750  mb  layer  and  is  characterized 
by  a  helical  type  circulation  about  a  horizontal  axis. 
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have  their  largest  amplitude  in  the  lowest  level.  Growth 
rates  for  wavenumbers  6  and  9  for  the  mean  winter  jet  and 
wavenumber  6  for  the  genesis  jet  are  reduced  by  less  than 
10%.  All  other  growth  rates  are  reduced  over  15%.  Growth 
rates  for  wavenumbers  greater  than  18  for  the  mean  winter 
jet  could  not  be  calculated  because  these  perturbati ons  were 
strongly  damped. 

The  most  unstable  waves  for  the  mean  winter  jet  in  the 
C  experiments  is  wavenumber  9  with  a  growth  rate  of  0.48 
day"^.  Therefore,  the  effect  of  friction  on  the  mean  winter 
jet  is  to  shift  the  most  unstable  wave  to  a  longer  wave¬ 
length  (wavenumber  12  in  set  A  versus  wavenumber  9  in  set  C 
experiment).  For  the  genesis  jet,  a  maximum  occurs  at  wave 
number  8  (1.04  day-1)  as  in  the  set  A  experiments  and  at 
wavenumber  36  (1.38  day-1).  The  growth  rate  spectrum  for 
the  genesis  jet  is  very  flat  for  wavenumbers  greater  than 
24,  which  suggests  that  the  helical  rolls  formed  in  the 
"mixed"  baroclinic  wave  tend  to  be  larger  when  frictional 
effects  are  included. 

The  growth  rates  for  the  modified  genesis  jet  with 
linearized  friction  (experiment  D)  are  shown  in  Fig.  3-21. 
Again,  the  shortest  wavelengths  a.e  most  highly  damped  with 
the  growth  rates  for  wavenumbers  greater  than  28  too  small 
to  calculate.  The  peak  in  the  growth  rate  spectrum  is  0.99 
day-1  at  wavenumber  8.  This  indicates  that  the  combination 
of  wave-zonal  flow  interaction  and  surface  friction  alter 


2.0 


experiments  (light  line). 


D.  EFFECT  OF  FRICTION 


An  experiment  (set  C)  similar  to  those  described  in  the 
previous  two  sections  is  performed  to  study  the  effect  of 
friction  on  the  linear  growth  rates.  These  C  experiments 
differ  from  the  set  A  experiments  in  that  the  surface 
friction  is  included.  The  zonal  mean  quantities  of  u,  v,  T 
and  it  are  again  restored  each  time  step. 

The  i ncorporati on  of  friction  in  the  C  experiments 
differs  from  the  form  described  in  Appendix  A.  Here,  the 
surface  stress  is 

t  =  P  C  q  | V | V  (3.3) 

where  CQ  =  l.OxlO-3.  It  is  assumed  that  the  surface  stress 
decreases  to  zero  linearly  with  pressure  through  the  lowest 
two  model  layers  (approximately  120  mb).  This  formulation 
of  the  surface  friction  is  used  to  provide  linear  friction 
over  the  entire  grid.  Although  this  form  is  proportional  to 
the  velocity  squared,  interactions  with  the  zonal  flow  are 
zero  since  the  zonal  flow  is  held  constant  in  time.  The 
i nteracti ons between  waves  is  small  since  the  perturbation 
i  s  smal 1 . 

The  major  effect  of  the  friction  is  to  decrease  the 
growth  rate  of  all  the  waves  (Fig.  3-20).  Also,  the 
shortest  waves  are  most  highly  damped,  which  is  consistent 
with  Gall  (1976b)  and  Gall  et  al .  (1979).  This  is  expected 
since  the  short  waves  have  been  shown  to  be  shallow  and  to 
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instability  creates  a  structure  different  from  the  observed 
structure  of  polar  lows.  Wave  G20  is  barotropi cal ly  damped 
since  HCK  <  0. 

The  results  for  the  genesis  jet  disagree  with  the 
findings  of  Sardie  and  Warner  (1983).  Using  a  3-layer 
quasi -geostrophi c  model,  they  found  the  most  unstable  wave¬ 
length  to  be  2700  km  (about  wavenumber  10)  with  a  growth 
rate  of  1.29  day-1.  They  did  not  find  large  growth  rates 
for  the  "mixed"  baroclinic  wave  (G40)  because  quasi- 
geostrophic  models  neglect  the  non-geostrophi c  effects 
necessary  for  Kel v i n-Hel mhol tz  instability.  Sardie  and 
Warner  (1983)  showed  that  inclusion  of  moist  processes 
shifted  the  most  unstable  wave  to  shorter  wavelengths  (600- 
2300  km  depending  on  the  heating  distribution)  and  concluded 
that  moist  baroclinic  processes  are  the  primary  cause  of 
Pacific  polar  lows. 

The  growth  rates  computed  in  this  study  for  the  modified 
genesis  jet  indicate  that  moist  processes  may  not  be  neces¬ 
sary  for  polar  lows  to  form.  Mullen  (1982)  documented  a 
case  of  deep  continental  polar  low  development  with  a  scale 
of  1800  km  in  a  domain  with  negligible  water  content.  The 
formation  of  a  polar  low  with  this  horizontal  scale  and 
without  moist  processes  is  consistent  with  the  findings  of 
this  study.  Some  other  process(es)  (possibly  terrain 
induced)  must  be  responsible  for  the  deep  structure  of  the 
continental  systems  observed  by  Mullen  (1979). 


Striking  similarities  exist  between  wavenumber  20  for 
the  modified  genesis  jet  (G20)  and  the  observations  of 
Pacific  polar  lows  by  Reed  (1979)  and  Mullen  (1979).  Reed 
(1979)  observed  that  the  polar  low  has  typical  wavelengths 
from  1000-1500  km.  Wave  G20  has  a  wavelength  of  1400  km. 
Mullen  (1979)  presented  composited  data  which  exhibited  R. 
values  of  1-2  at  the  position  of  the  polar  low  from  the 
surface  to  850  mb.  This  has  also  been  found  to  be  the  case 
for  G20  (cf,  Fig.  3-16).  Mullen  (1979)  also  showed  that 
the  polar  low  formed  500-700  km  north  of  the  jet  while  G20 
formed  about  1000  km  north  of  the  jet.  Differences  exist 
between  the  vertical  structure  of  observed  polar  lows  and 
G20.  Reed  (1979)  documented  a  Pacific  polar  low  which 
exhibited  geopotential  perturbations  at  300  mb  as  large 
as  those  observed  at  the  surface.  G20  has  been  found  to 
exhibit  maximum  amplitude  at  the  surface  with  little 
evidence  of  the  wave  above  500  mb.  It  has  been  shown  by 
Gall  (1976c)  that  latent  heat  release  can  serve  to  deepen 
systems.  This  process  will  be  examined  in  Chapter  IV. 

It  has  been  demonstrated  that  G20  is  formed  by  conven¬ 
tional  baroclinic  instability  which  supports  the  hypotheses 
of  Reed  (1979)  and  Mullen  (1979)  that  Pacific  polar  lows  are 
barocl i ni cal  1 y  unstable  disturbances.  Energy  conversion 
processes  suggest  Kel vi n-Helmhol tz  instability  is  a  con¬ 
tributing  factor,  but  it  is  about  four  times  smaller  in 
magnitude  than  CE.  In  addition,  the  Kel vi n-Helmhol tz 
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Fig.  3-19.  Surface  pressure  at  84h  for  wavenumber  20 
in  modified  genesis  jet.  Contour  interval  is 
4  mb.  Figure  represents  a  subsection  of  the 
forecast  domain.  Southern  boundary  of  figure 
is  at  the  center  of  the  forecast  domain. 

100  km  of  right  side  of  domain  is  the  same 
as  leftmost  100  km  due  to  periodicity. 
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3-18.  Similar  to  Fig.  3-10  except  for  wavenumber  20 
in  modified  genesis  jet.  Contour  interval  is 
2  ubar/s  for  w  and  10  K  for  T.  100  km  of  right 
side  of  domain  is  the  same  as  the  leftmost 
100  km  due  to  periodicity. 
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Similar  to  Fig.  3-7  except  for  wavenumber  20 
in  modified  genesis  jet. 


The  deepening  rate  at  any  latitude  <j> ,  is  adjusted  by  multi¬ 
plying  by  sin  60°/sin  <t> .  Here  G20D  exhibits  a  13  mb  fall 
in  24  h  near  52N,  which  corresponds  to  a  deepening  rate  of 
0.6  B.  The  maximum  deepening  rate  for  G20A  is  0.22  B. 

Mullen  (1983)  showed  that  Pacific  polar  lows  can  achieve 
deepening  rates  of  about  1.9  B,  although  this  seems  to  occur 
only  when  the  polar  low  interacts  with  the  main  frontal  band 
equatorward  of  the  deepening  low.  For  a  polar  low  which 
did  not  interact  with  the  main  polar  front,  Mullen  found  a 
deepening  rate  of  0.85  B.  This  corresponds  well  with  the 
deepening  rate  computed  here  for  G20D.  Observed  deepening 
rates,  such  as  those  computed  by  Mullen,  should  be  expected 
to  be  larger  than  those  computed  in  this  study.  Pacific 
polar  lows  generally  deepen  as  they  move  into  a  large  area 
of  low  pressure  often  centered  in  the  Gulf  of  Alaska. 
Therefore,  some  of  the  pressure  fall  must  be  attributed 
to  the  new  environment  of  the  storm,  rather  than  actual 
intensification  of  the  individual  cyclone.  This  "change  of 
environment"  did  not  exist  in  G20D  and  therefore,  deepening 
rates  computed  using  real  data  could  very  well  be  larger 
than  those  computed  here. 

The  sea-level  pressure  fields  at  60,  72,  84  and  96  h  for 
G20A  and  G20D  are  shown  in  Figs.  4-2  and  4-3,  respectively. 
During  intensification,  the  low  center  in  G20D  (hereafter 
G20D)  moves  due  east  while  G20A  moves  toward  the  east- 
northeast.  At  96  h,  the  center  in  G20A  is  about  200  km 
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Sea-level  pressure  fields  for  G20A  at  (a) 

60  h  (b)  72  h,  (c)  84  h  and  (d)  96  h.  Figures 
represent  a  subsection  of  forecast  domain. 
Southern  boundary  of  figure  is  located  200  km 
north  of  center  of  forecast  domain.  100  km  of 
right  side  of  domain  is  the  same  as  the  leftmost 
100  km  due  to  periodicity. 


north  of  the  correspondi ng  position  of  G20D.  Since  diabatic 
processes  are  included  in  G20D,  it  is  suspected  that  the 
sea-surface  temperature  (Fig.  4-4)  plays  a  role  in  deter¬ 
mining  the  position  of  the  low.  Sanders  and  Gyakum  (1980) 
demonstrated  that  strong  oceanic  development  tends  to  occur 
near  the  maximum  gradient  of  the  sea-surface  temperature. 
Here,  it  is  clear  that  during  the  developing  stages,  the 
center  in  G20D  remains  in  the  region  of  maximum  gradient  of 
sea-surface  temperature  (near  the  6°C  isotherm),  while  the 
center  in  G20A  moves  away  from  it. 

The  low  in  G20A  is  observed  to  form  about  900-1000  km 
north  of  the  jet  axis  while  the  low  in  G20D  forms  about 
750  km  north  of  the  jet  axis.  This  is  consistent  with  the 
composite  results  of  Mullen  (1979),  which  indicated  that 
the  polar  low  forms  500-700  km  north  of  the  jet.  Therefore, 
realistic  deepening  rates  and  more  realistic  storm  positions 
are  produced  when  diabatic  effects  are  included. 

Both  G20A  and  G20D  exhibit  a  band  of  positively  corre¬ 
lated  u'  and  oo'  similar  to  that  in  G40  (cf,  Fig.  3-13) 
parallel  to  the  trough  extending  from  the  low  center.  In 
this  region,  the  <  1,  which  suggests  that  Kelv in- 
Helmholtz  instability  may  be  significant  near  the  trough. 

East-west  cross-sections  of  height  perturbat i ons  at  84  h 
are  given  in  Fig.  4-5  for  G20A  and  G20D.  These  figures  are 
representative  of  the  developing  stages  of  the  disturbance. 
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G20A  exhibits  a  maximum  deviation  at  the  surface  of  about 
16  m.  The  wave  is  damped  strongly  in  the  vertical  so  that 
the  amplitude  at  800  mb  is  only  50%  of  the  surface  value. 
620D  is  noticeably  stronger  with  a  surface  maximum  amplitude 
of  110  m.  The  amplitude  does  not  decrease  as  rapidly  in  the 
vertical  as  in  G20A.  The  amplitude  at  500  mb  is  about  one- 
half  the  surface  amplitude. 

Since  G20D  extends  through  a  greater  depth  than  G  2  0  A ,  it 
should  be  expected  to  have  a  phase  speed  greater  than  that 
of  G20A.  Phase  speeds  computed  from  the  surface  pressure 
centers  in  Figs.  4-2  and  4-3  are  12.4  m  s-1  for  G20D  and 
8.2  m  s-1  for  G20A.  Some  of  the  increase  in  phase  speed  of 
G20D  must  also  be  due  to  it  being  imbedded  in  stronger  flow 
along  the  more  southern  track. 

East-west  cross-sections  of  vertical  motion  and  temper¬ 
ature  at  84  h  are  shown  in  Fig.  4-6.  The  reason  for  the 
increased  vertical  scale  in  G20D  becomes  obvious  upon 
inspection  of  the  u  and  T  fields.  For  G20A,  the  warmest 
air  extends  upward  from  the  surface  to  about  750  mb  with 
little  or  no  tilt.  The  maximum  vertical  motion  occurs  just 
east  of  the  warmest  air,  which  favors  generation  of  KE. 
Above  750  mb,  the  warmest  air  tilts  to  the  west  with  height 
away  from  the  maximum  vertical  velocity  similar  to  the 
pattern  exhibited  by  wavenumber  20  in  the  linear  experiments 
in  Chapter  III.  In  addition,  the  temperature  perturbation 
decreases  above  750  mb.  These  factors  are  consistent  with  a 


reduced  generation  or  even  a  destruction  of  KE  above  750  mb. 
This  will  be  discussed  further  in  the  next  section.  The 
structure  of  G20A  strongly  resembles  that  of  G20  in  the 
linear  experiments  in  Chapter  III.  This  is  not  surprising 
since  no  other  waves  are  included  in  these  studies,  and 
since  the  disturbance  is  relatively  weak,  wave-zonal  flow 
interaction  is  small.  In  contrast  to  the  vertical  structure 
of  G20A,  the  low  in  G20D  exhibits  a  strong  positive  tempera¬ 
ture  maximum  which  extends  from  the  surface  to  450  mb  with 
little  or  no  tilt.  This  results  in  a  generation  of  KE  over 
a  much  deeper  layer  than  in  G20A. 

One  reason  for  the  increased  depth  of  the  temperature 
perturbation  in  G20D  can  be  seen  by  examining  Fig.  4-7, 
which  shows  the  amount  of  total  ,  convective  and  stable 
precipitation  that  accumulated  during  the  12  h  period  ending 
at  84  h.  A  maximum  of  precipitation  (about  22  mm)  is  found 
at  and  just  ahead  of  the  path  of  the  low  center  with  a 
secondary  maximum  (about  15  mm)  south  of  the  low  along  the 
front.  About  70%  of  the  precipitation  occurring  near  the 
low  center  and  about  85%  in  the  southern  maximum  is  convec¬ 
tive.  Based  on  cross-sections  of  9g  and  relative  humidity 
(not  shown),  conditional  instability  extended  up  to  about 
500  mb  while  stratiform  clouds  (relative  humidity  100%) 
extended  up  to  about  350  mb  in  the  warm  sector  ahead  of  the 
low.  The  release  of  latent  heat  over  this  depth  is  consist¬ 
ent  with  the  warming  observed  in  G200  ahead  of  the  low.  The 
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Fig.  4-7.  1 2 h  accumulation  of  (a)  total,  (b)  convective 

and  (c)  stable  precipitation  in  a  fixed-frame 
reference  ending  at  84  h  for  G200.  Contour 
interval  is  5  mm.  Grid  placement  as  in 

Fig.  4-2  . 
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diabatic  heating  effectively  reduces  the  amount  of  adiabatic 
cooling  in  the  region  of  ascent. 

To  separate  the  effects  of  latent  heat  release  by 
precipitation  from  the  PBL  diabatic  processes,  a  simulation 
(G20E)  was  made  which  was  similar  to  G20D  except  that  latent 
heat  release  by  precipitation  was  not  allowed.  G20E  attains 
a  minimum  central  pressure  of  993  mb  (Fig.  4-8),  which  is 
10  mb  lower  than  the  adiabatic  simulation,  but  9  mb  higher 
than  the  fully  diabatic  simulation.  Obviously,  latent  heat 
release  by  precipitation  is  not  acting  alone  in  the  intensi¬ 
fication  of  G20D.  Significant  contributions  to  the 
intensification  of  the  polar  low  are  made  by  three  other 
processes.  First,  the  sensible  heating  due  to  the  large 
air-sea  temperature  difference  insures  a  strong  meridional 
low-level  temperature  gradient  (i.  e.,  baroclinity)  over  a 
longer  time  period  than  found  in  the  adiabatic  simulation. 
Second,  the  sensible  heating  within  the  PBL  can  be  realized 
in  the  mid-levels  almost  immediately  through  dry  convective 
adjustment.  This  is  because  the  lapse  rate  is  nearly  dry 
adiabatic  in  the  genesis  jet  from  the  surface  to  near  750  mb 
in  the  vicinity  of  the  polar  low.  Third,  the  large,  initial 
air-sea  temperature  difference  acts  to  maintain  a  nearly  dry 
adiabatic  lapse  rate  in  the  lowest  troposphere  ( i .  e . ,  low 
R i )  which  should  increase  the  growth  rates  of  the  shorter 
waves  (Staley  and  Gall,  1977).  To  show  the  influence  of  the 
surface  heat  flux,  G20E  was  rerun  without  a  perturbation  in 
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e  initial  fields.  The  top  of  the  low  level  dry  adiabatic 
yer  near  the  latitude  of  the  polar  low  rose  from  near 
0  mb  initially  to  near  650  mb  by  96  h.  The  effect  of 
ese  two  processes  on  the  energetics  of  G20D  will  be 
scussed  further  in  the  following  section. 

The  effect  of  the  diabatic  processes  on  G20  can  be  seen 
Figs.  4-9  and  4-10,  which  show  the  height  and  temperature 
elds  at  84  h  for  850,  700  and  500  mb  for  G20A  and  G20D, 
spectively.  The  amplitude  of  G20D  is  stronger  at  all 
vels.  In  addition,  the  warmest  (coldest)  air  in  G20A 
mains  just  to  the  east  of  the  trough  (ridge)  position  at 
1  levels.  The  amplitude  of  the  thermal  ridge  is  larger 
d  farther  to  the  east  of  the  trough  axis  in  G20D,  espe- 
al ly  at  700  and  500  mb. 

The  surface  sensible  and  latent  heat  fluxes  for  G20D  at 
h  are  shown  in  Fig.  4-11.  The  maximum  for  both  is  about 
0  W  m-2  and  is  found  to  the  southwest  of  the  low  center 
hind  the  cold  front.  This  occurs  because  this  quadrant 
hibits  relatively  low  temperatures  and  the  strongest  low- 
vel  winds.  Although  downward  fluxes  are  typically 
served  in  the  warm  sector  ahead  of  a  low,  G20D  exhibits  a 
condary  maximum  of  upward  surface  fluxes  roughly  one-half 
e  maximum  values  in  the  cold  air.  There  are  two  reasons 
r  the  upward  fluxes  in  the  warm  air.  First,  the  initial 
r-sea  temperature  difference  in  G200  was  at  least  5°C  in 
e  vicinity  of  the  developing  low,  a  condition  that  results 
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Fig.  4-9.  Fields  of  height  and  temperature  at  (a)  850  mb, 
(b)  700  mb  and  (c)  500  mb  at  84  h  for  G20A. 
Contour  interval  for  height  is  30  m  in  (a) 
and  (b)  and  60  m  in  (c).  Contour  interval 
for  temperature  is  5°C.  Grid  placement  as  in 
Fig.  4-2. 
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Fig.  4-10.  Similar  to  Fig.  4-9  except  for  G20D 


heat  flux  at  84  h  for  G20D.  Contour  interval  is 
100  W  m-2.  Grid  placement  as  in  Fig.  4-2. 


in  an  upward  surface  heat  flux  in  all  quadrants  around  the 
low.  Second,  the  low-level  horizontal  temperature  perturba¬ 
tion  associated  with  G200  (Fig.  4 - 1 0 a )  did  not  become  strong 
enough  (at  least  through  84  h)  to  reverse  the  sign  of  the 
initial  air-sea  temperature  difference. 

B.  ENERGETICS 

The  equations  of  the  energy  budgets  used  in  this  section 
are  given  in  Appendix  B.  However,  the  terms  involving 
changes  in  the  surface  pressure  (i.  e.,  PAE,  PAZ,  PKE,  PKZ) 
are  small  relative  to  the  other  terms,  and  therefore  are  not 
included  in  this  discussion.  The  budgets  are  not  perfectly 
balanced  due  to  truncation,  interpolation,  round-off  errors, 
etc.  For  the  adiabatic  runs,  the  residual  is  generally  an 
order  of  magnitude  less  than  the  leading  term(s),  and  the 
residual  for  the  diabatic  runs  is  about  5  times  less  than 
the  leading  term(s). 

KE  budgets,  integrated  over  the  forecast  domain,  are 
presented  in  6  h  intervals  in  Tables  4-1  and  4-2  for  G20A 
and  G20D,  respectively.  Steady  increases  of  KE  are  noted  in 
both  runs  from  the  beginning  of  the  forecast.  The  maximum 
KE  attained  in  G20A  (7  x  1Q3  J  m~2)  is  an  order  of  magnitude 
less  than  the  maximum  reached  in  G20D  (8  x  10^  J  m-2).  In 
addition,  the  maximum  value  occurred  from  114-132  h  for  G20A 
and  at  90  h  for  G200  ,  which  shows  the  importance  of  the 
diabatic  processes  in  increasing  both  the  deepening  rate  as 
well  as  the  intensity  of  the  system. 
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Table  4-1.  KE  budgets  for  G20A.  KE  in  units  of  1 D 2  Jm 


A 1  1 

otHer  budget 

terms 

in  units 

of  TO-' 

Wm-2  . 

Time  (hr) 

KE 

aKE/at 

CE 

HCK 

VCK 

DE 

48 

2.4 

.05 

.04 

.00 

.03 

-  .03 

54 

3.5 

.07 

.07 

-.01 

.04 

-  .04 

60 

5.3 

.10 

.10 

-.01 

.06 

-  .06 

66 

8.0 

.15 

.16 

-.02 

.08 

-  .10 

72 

11.7 

.21 

.23 

-.04 

.12 

-  .14 

78 

17.1 

.29 

.34 

-.07 

.18 

-  .21 

84 

24.0 

.36 

.46 

-.11 

.25 

-  .31 

90 

32.8 

.44 

.65 

-.18 

.34 

-  .46 

96 

42.8 

.46 

.82 

-.27 

.44 

-  .65 

102 

52.5 

.42 

.98 

-.38 

.53 

-  .87 

108 

61.0 

.33 

1.09 

-.48 

.62 

-1.-08 

114 

66.9 

.21 

1.12 

-.55 

.69 

-1.25 

1  20 

69.9 

.08 

1.06 

-.57 

.77 

-1.34 

126 

70.3 

.00 

.93 

-.50 

.82 

-1.36 

132 

70.1 

.08 

.86 

-.34 

.87 

-1.34 

138 

73.9 

.29 

1.01 

-.24 

.90 

-1.36 

100 


t 

Table 

4-2.  S i m i lar  to 

Table  4-1 

except 

for  G20D . 

t 

Time  (hr) 

KE 

a  KE/at 

CE 

HCK 

VCK 

OE 

vK; 

48 

63.7 

1.55 

1.0 

.1 

1.5 

.3 

\ 

54 

94.1 

1.48 

1.2 

-.5 

1.8 

-.5 

60 

127.6 

2.13 

1.6 

-.7 

2.7 

-1.0 

66 

186.2 

3.34 

3.7 

-1.6 

4.5 

-1.5 

72 

272.0 

5.22 

4.6 

-1.3 

5.1 

-1.9 

78 

411.8 

8.46 

10.8 

-1.6 

3.8 

-3.8 

84 

637.6 

9.15 

16.6 

-4.4 

6.3 

-7.2 

90 

807.2 

1.78 

15.6 

-13.5 

12.4 

-9.3 

96 

714.6 

-6.45 

4.3 

-13.0 

10.5 

-7.8 

v_  V- 

102 

528.4 

-7.20 

1.5 

-6.1 

4.9 

-5.4 

■ 

108 

394.8 

-3.58 

1.0 

-2.5 

4.1 

-3.6 

*.'*’-* 

114 

373.6 

-1.25 

3.6 

-3.6 

3.1 

-3.9 

120 

341.0 

-.84 

5.0 

-4.2 

2.3 

-3.0 

-  ..4 

126 

337.5 

1.27 

3.7 

-1.4 

2.2 

-2.3 

132 

396.0 

2.55 

3.4 

-1.4 

4.0 

-2.9 

>v: 

138 

447.8 

.26 

2.2 

-3.5 

5.3 

-4.1 
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Time-height  cross-sections  of  KE  for  G20A  and  G20D  are 
given  in  Fig.  4-12.  For  G20A,  the  maximum  KE  occurs  at  950 
mb  and  decreases  rapidly  with  height.  For  G20D,  the  maximum 
KE  is  found  at  90  h,  also  near  950  mb,  but  a  secondary  maxi¬ 
mum  of  about  2/3  of  the  950  mb  value  exists  near  450  mb. 

This  is  a  reflection  of  the  increased  depth  of  this  system, 
as  noted  in  the  previous  section. 

During  the  growth  of  G20A  and  G20D,  CE  is  the  largest 
term,  which  indicates  that  baroclinic  conversion  is  the 
principal  mechanism  in  forming  the  1400  km  wave.  The  values 
of  CE  are  more  than  an  order  of  magnitude  larger  for  G20D 
than  for  G20A.  Time-height  cross-sections  of  CE  are  shown 
in  Fig.  4-13.  For  G20A,  CE  has  a  maximum  at  120  h  near 
37  5  mb  with  negative  values  above  700  mb  as  the  storm  is 
developing.  This  is  consistent  with  the  structure  of 
G20A  found  in  the  previous  section.  It  was  noted  that  a 
thermally  direct  circulation  was  evident  below  about  750  mb, 
which  is  consistent  with  the  shallow  nature  of  CE  in  Fig. 
4-13.  For  G20D,  CE  exhibits  a  maximum  at  34  h  near  700  mb 
that  is  nearly  an  order  of  magnitude  larger  than  the  G20A 
maximum.  In  addition,  positive  values  of  CE  extend  to 
nearly  400  mb.  It  was  noted  in  the  previous  section  that 
the  effect  of  diabatic  heating  was  to  shift  the  warmest  air 
above  750  mb  to  the  region  of  the  strongest  upward  motion, 
thereby  enabling  the  generation  of  KE,  through  the  term  CE, 
to  strengthen  the  circulation  in  the  middle  to  upper  levels. 


ght  solid  line  is  positive  cont 
gative  contour. 


The  horizontal  shear,  or  the  barotropic  term,  HCK,  is 
about  one  order  of  magnitude  less  than  CE  for  both  G20A  and 
G20D  during  the  early  stages  of  development.  As  the  wave 
grows,  HCK  is  negative  and  increases  in  magnitude.  In  G20A, 
HCK  reaches  a  value  of  about  one-half  that  of  CE  at  120  h 
while  in  G20D,  HCK  becomes  the  largest  term  at  96  h.  There¬ 
fore,  barotropic  effects  are  important  in  halting  the  growth 
of  the  polar  low,  especially  when  diabatic  effects  are 
included.  Time-height  cross-sections  of  HCK  are  given  in 
Fig.  4-14.  As  with  CE,  G20A  exhibits  a  maximum  of  HCK  near 
900  mb  while  G20D  exhibits  a  maximum  at  850  mb  and  near 
500  mb. 

The  Kelvin-Helinholtz  term,  VCK,  is  positive  at  all  times 
for  both  forecasts.  It  was  shown  in  Chapter  III  that  VCK 
involves  the  correlation  between  u'  and  w1  multiplied  by  the 
vertical  wind  shear.  During  the  development  of  the  polar 
low  in  these  experiments,  strong  positive  correlations  were 
observed  between  u'  and  go'.  Generally,  this  implies 
easterly,  upward  motion  ahead  of  the  trough  and  westerly, 
downward  motion  behind  the  trough.  Coupled  with  the  large 
vertical  wind  shear,  large,  positive  values  of  VCK  result. 

It  was  shown  in  the  previous  section  that  the  sensible  heat 
flux  in  G20D  reduces  the  stability  in  the  lowest  levels, 

which  results  in  a  decrease  of  Ri»  which  is  consistent  with 
the  large  values  of  VCK.  However,  VCK  is  ruled  out  as  the 
term  which  is  responsible  for  the  rapid  deepening  of  the 
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polar  low.  In  fact,  VCK  may  inhibit  polar  low  formation. 

As  described  in  Chapter  III,  positive  values  of  VCK  lead  to 
helical  circulations  in  the  vertical  plane  which  are  quite 
different  from  the  observed  structure  of  polar  lows.  The 
presence  of  this  helical  circulation  will  tend  to  decrease 
the  vertical  wind  shear  (i.e.,  baroclinity)  and  result  in 
a  reduced  growth  rate  of  the  polar  low  through  CE.  Also, 
during  the  rapid  growth  of  G200  from  78-90  h,  CE  is  the 
dominant  term.  For  example,  VCK  is  three  times  less  than 
CE  at  78  n.  When  weaker  development  was  occurring  prior 
to  78  h,  VCK  was  larger  than  CE  in  G200.  In  G20A,  VCK  was 
generally  a  factor  of  two  smaller  than  CE  at  all  times. 
Time-height  cross-sections  of  VCK  are  shown  in  Fig.  4-15. 
Again,  the  increased  depth  of  the  disturbance  is  reflected 
by  a  double  maximum  in  VCK  at  850  mb  and  500  mb  in  G20D  and 
only  a  low-level  maximum  of  VCK  near  850  mb  for  G20A. 

The  dissipation  of  KE  by  friction,  DE,  is  negative  at 
all  t’ines,  which  indicates  a  sink  of  KE.  For  G20A,  DE 
slow  ^  s  e  s  in  magnitude  during  the  forecast  so  that 

at  24  h,  DE  is  one  half  that  of  CE.  At  114  h,  DE  is  larger 
than  CE,  which  shows  the  strong  influence  of  friction  in 
damping  this  disturbance  at  all  times.  In  fact,  DE  is  the 
largest  term  in  the  G20A  budget  from  114-138  h.  This  is 
expected  due  to  the  shallow  nature  of  this  disturbance. 

For  G20D,  DE  plays  a  reduced,  but  still  significant  role. 
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The  results  of  experiments  G20B,  G20C  and  G20D 
illustrate  the  type  of  environment  that  is  favorable  for 
significant  polar  low  development.  As  the  polar  low 
develops,  relatively  cold  air  is  advected  to  the  west  of 
the  storm  in  the  lower  levels,  while  warmer  air  is  advected 
ahead  of  the  storm.  This  temperature  advection  drives  the 
PBL  toward  more  unstable  conditions  behind  the  storm  and 
toward  more  stable  conditions  ahead  of  the  storm.  For  an 
initially  neutral  PBL  (as  in  G20B),  this  implies  an  unstable 
PBL  in  the  cold  air  and  a  stable  PBL  in  the  warm  air  during 
the  growth  stage.  Therefore,  surface  fluxes  in  the  warm 
sector  tend  to  be  relatively  small.  This  implies  that 
relatively  little  or  no  additional  moisture  will  be  added  to 
the  system  which  may  be  used  to  sustain  additional  growth. 

In  addition,  surface  friction  would  also  be  relatively  small 
in  the  warm  sector  which  implies  smaller  vertical  motion  and 
smaller  moisture  convergence,  and  therefore  a  reduced  amount 
of  precipitation.  In  the  unstable  PBL  behind  the  surface 
low,  relatively  large  surface  fluxes  exist.  However, 
heating  in  this  area  will  destroy  AE.  More  favorable  condi¬ 
tions  for  development  exist  for  an  initially  unstable  PBL  in 
all  quadrants  around  the  developing  disturbance  ( G 2 0 D ) .  In 
this  case,  large  amounts  of  heat  and  moisture  are  added  to 
the  near-surface  layers,  which  is  a  favorable  condition  for 
convection.  In  addition,  surface  friction  is  increased, 


1010 


*o 

V 


o 

-a 


CO 

o 

C\J 

CJ3 


(awl  aanssaad 


ca 

co 

CD 


The  minimum  central  pressures  versus  time  for  all  four 
G20  simulations  are  given  in  Fig.  4-20.  The  least  intense 
storm  was  produced  by  the  adiabatic  simulation,  G20A,  which 
attained  a  minimum  pressure  of  1003  mb.  Although  this 
forecast  did  not  include  diabatic  effects  explicitly,  the 
initial  air-sea  temperature  difference  did  play  a  role  in 
the  evolution  of  the  storm  since  the  surface  friction  is  a 
function  of  the  stability  of  the  PBL  (Eqs.  A-10  through 
A-20).  Initially,  the  PBL  was  unstable  due  to  the  large 
air-sea  temperature  difference  imposed  and  remained  unstable 
throughout  the  forecast.  This  is  somewhat  unrealistic 
because  the  surface  heat  flux  would  act  to  drive  the  PBL 
toward  a  neutral  stability.  Therefore,  G20A  was  rerun  using 
an  initial  air-sea  temperature  difference  of  zero  to  deter¬ 
mine  the  effect  of  PBL  stabilization  on  G20A.  Although  no 
difference  in  storm  structure  or  speed  was  noted,  the 
minimum  central  pressure  did  decrease  2  mb  more  than  in 
G20A. 

The  simulation  of  G20D  produces  the  most  intense  storm 
and  the  storm  with  the  greatest  vertical  depth.  The  rapid 
deepening  of  the  system  is  crucially  tied  to  the  latent  heat 
release.  From  72-84  h,  when  G20C  and  G20D  exhibit  the  most 
significant  deepening,  there  is  15-22  mm  of  precipitation 
near  and  ahead  of  the  storm  center  in  G20D.  During  the 
same  time,  G20C  has  5-7  mm  and  G20B  has  less  than  3  mm  of 
precipitation. 
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tends  to  be  too  shallow.  Inclusion  of  diabatic  processes 
leads  to  a  significant  increase  in  intensity  (i.e.,  19  mb 
lower  central  pressure)  and  a  deeper,  more  realistic  verti¬ 
cal  structure. 

C.  EFFECT  OF  AIR-SEA  TEMPERATURE  DIFFERENCE' 

In  the  previous  two  sections,  large  differences  were 
observed  between  the  adiabatic  and  diabatic  experiments. 

The  air-sea  temperature  differences  given  by  Mullen  (1979) 
for  the  genesis  case  (cf.  Fig.  4-4)  were  used  in  the  initial 
fields.  The  results  of  two  additional  simulations  of  G20 
are  presented  here  to  illustrate  the  effect  of  the  air-sea 
temperature  difference.  Details  of  all  four  G20  simulations 
are  given  in  Table  4-5.  A  zero  initial  air-sea  temperature 
difference  is  specified  in  G20B,  while  in  G20C  the  air-sea 
temperature  difference  is  one-half  of  that  used  in  G20D. 

All  diabatic  processes  used  in  G20D  are  also  used  in  G20B 
and  G20C. 

Table  4-5.  Physical  processes  and  initial  air-sea  tempera¬ 
ture  differences  for  polar  low  simulations. 

Initial  Air-Sea 


Experi ment 

Friction 

Diabatic 

Processes 

Temperature 
Difference  ( K ) 

G20  A 

Yes 

No 

Mullen  (1979) 
Genesis  Case 

G20B 

Yes 

Yes 

0.0 

G20C 

Yes 

Yes 

One-half  value 
used  in  G20A 

G20D 

Yes 

Yes 

Same  as  in  G20A 

I 
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» 


9 
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The  most  significant  differences  between  the  forecasts  is 
the  large  increase  in  VCA  near  600  mb  in  G20D,  which  again 
reflects  the  increased  circulation  in  the  middle  layers. 

The  results  of  these  energy  budgets  illustrate  how  the 
diabatic  effects  lead  to  an  increase  in  storm  intensity. 
During  the  rapid  intensification  of  G20D  from  73-90  h,  the 
sum  of  HCA  and  VCA  is  6-10  times  larger  than  GE.  This 
suggests  that  much  more  AE  is  formed  by  eddy  heat  transport 
processes  than  generated  directly  by  diabatic  effects. 

This  is  supported  by  the  fact  that  experiment  G20E,  which 
neglected  latent  heat  release  by  precipitation  resulted  in 
a  10  mb  decrease  in  minimum  central  pressure  from  the  adia¬ 
batic  simulation.  Therefore,  it  is  concluded  that  diabatic 
processes  in  the  Pacific  polar  low  primarily  serve  to  alter 
the  overall  temperature  structure  of  the  disturbance  such 
that  eddy  heat  transport,  and  therefore,  AE  values  are 
increased.  The  rapid  storm  growth  is  due  to  baroclinic 
conversion  of  this  AE  to  KE. 

The  results  presented  here  agree  with  the  speculation  of 
Sardie  and  Warner  (1983)  that  moist  baroclinity  is  the  cause 
of  Pacific  polar  lows.  The  linear  studies  in  Chapter  III 
demonstrated  that  wavenumber  20  exhibits  growth  rates  as 
large  as  wavenumber  8  after  a  period  of  modification  to  the 
genesis  jet.  Although  wavenumber  20  forms  in  the  nonlinear, 
adiabatic  model  with  friction,  the  intensity  of  the  polar 
low  is  much  weaker  than  observed.  Furthermore,  the  wave 


-19.  Time-height  cross-sections  of 
G20D.  Contour  interval  is  5  x 
30  x  10" 3  W  m -2  in  (b)  Format 
Fig.  4-13. 
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Fig.  4-17.  Time-height  cross- secti on  of  GE  for  G20D. 

Contour  interval  is  10  x  1 0 - 3  W  m-2.  Format 
of  contours  as  in  Fig.  4-13. 
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GE  for  G200  is  shown  in  Fig.  4-17.  Maximum  values  of  GE  are 
found  near  800  mb  with  positive  values  extending  from  near 
950  mb  to  400  mb.  Small  negative  values  are  observed  near 
the  surface  due  to  the  largest  surface  fluxes  occurring  in 
the  cold  air  behind  the  low  (Fig.  4-11). 

The  horizontal  eddy  heat  transport  term,  HCA,  is 
positive  during  the  development  of  G20A  and  G200,  which 
indicates  a  generation  of  AE.  Maximum  values  of  HCA  for 
G20A  are  found  near  950  mb  with  smaller,  negative  values 
above  700  rnb  (Fig.  4-18).  For  the  diabatic  run,  HCA  is 
larger  in  magnitude  and  the  positive  values  extend  up  to 
450  mb.  This  illustrates  that  the  diabatic  processes  have 
resulted  in  an  increased  circulation  above  750  mb  such  that 
horizontal  eddy  heat  transport  can  generate  significant  AE 
at  these  levels.  This  is  expected  due  to  the  height  and 
temperature  patterns  exhibited  by  G20A  and  G20D  (Figs.  4-9 
and  4- 10 ) . 

The  vertical  component  of  CA,  VCA,  is  also  positive 
during  the  development  of  the  polar  low  in  both  experiments. 
Again,  G20D  exhibits  values  over  one  order  of  magnitude 
larger  than  G20A.  Although  the  term  containing  E'oj'ct  favors 
destruction  o-f  AE,  the  terms  containing  u>'E'  3T/5P  and 
(jj'T'aE/ap  are  positive  and  larger  in  magnitude,  and  result 
in  an  overall  gain  in  AE  through  VCA.  Time-height  cross- 
sections  of  VCA  for  G20A  and  G20D  are  given  in  Fig.  4-19. 
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Table  4-3. 

AE  budgets 

for  G20A . 

AE  in 

units  of  10 

2  Jm'2. 

-V J 

A1 1  other 

budget  terms 

in  units  of  1  0  " 1 

Wm”2 . 

‘j 

I  Time  (hr) 

AE 

3  A  E  /  3  t 

CE 

HCA 

VCA 

.  ■-  « 
:  •  ", 

l 

48 

1.3 

.03 

.04 

.03 

.04 

54 

1.9 

.04 

.07 

.05 

.06 

|  60 

2.9 

.05 

.10 

.08 

.09 

» 

66 

4.2 

.08 

.15 

.11 

.13 

72 

6.1 

.10 

.23 

.16 

.19 

78 

8.7 

.14 

.34 

.23 

.27 

*  * 

• 

84 

12.0 

.17 

.46 

.31 

.36 

90 

15.9 

.19 

.65 

.41 

.49 

!  96 

20.0 

.18 

.82: 

.50 

.58 

;  102 

23.7 

.14 

.98 

.57 

.65 

108 

26.1 

.06 

1.09 

.61 

.67 

I  114 

26.5 

-.02 

1.12 

.61 

.60 

120 

25.2 

-.05 

1.06 

.60 

.52 

126 

24.2 

.01 

.93 

.63 

.43 

|  132 

25.8 

.14 

.86 

.75 

.40 

• 

138 

30.3 

.26 

1.01 

.95 

.47 

*  . 

t 

-  4 
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Although  DE  attains  a  maximum  at  90  h  in  G20D  when  KE  is  a 
maximum,  0E  is  still  30%  smaller  than  H ' " .  Therefore,  when 
diabatic  processes  are  included,  barotropic  damping  is 
dominant  in  halting  wave  growth.  This  is  due  to  the 
increased  depth  of  G20D.  Frictional  effects  are  relatively 
small  in  G200  since  they  are  felt  only  in  the  lowest  levels 
of  the  model.  However,  it  was  shown  earlier  (Fig.  4-14) 
that  HCK  exhibits  large  values  near  850  mb  and  near  500  mb 
and  is,  therefore,  more  efficient  in  reducing  overall  wave 
growth . 

AE  budgets,  integrated  over  the  forecast  domain,  are 
presented  in  6  h  intervals  in  Tables  4-3  and  4-4  for  G20A 
and  G20D,  respectively.  G20D  exhibits  more  AE  than  G20A 
at  all  times.  The  maximum  AE  for  G20A  (2.6  x  1 0 3  j  m-2) 
is  found  at  114  h,  and  the  maximum  AE  for  G20D  (1.3  x 
10^  J  m-2)  is  found  at  90  h.  Therefore,  as  evidenced  with 
KE,  AE  develops  faster  and  reaches  larger  values  with  the 
inclusion  of  diabatic  processes.  Time-height  cross-sections 
of  AE  are  given  in  Fig.  4-16.  The  major  difference  between 
the  forecasts  is  that  G20D  exhibits  a  strong  maximum  of  AE 
near  300  mb.  This  feature  is  absent  in  the  G20A  forecast. 

The  generation  of  AE  by  diabatic  processes,  GE,  is 
positive  at  all  times  in  G20D  as  a  result  of  the  precipita¬ 
tion  found  in  advance  of  the  storm  center  (cf,  Fig.  4-8). 

The  largest  values  of  GE  occur  from  78-90  h  during  the  rapid 
intensification  of  G200.  The  time-height  cross-section  of 
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Contour  interval  is  b  x  iu-J  w  m-^  in  t  a ; 
10-3  w  m-2  in  (b).  Format  of  contours  as  in 
4-13. 


which  suggests  stronger  vertical  motion  and  therefore,  an 
increase  in  precipitation.  Although  AE  is  still  destroyed 
in  the  cold  air  behind  the  low,  this  is  small  in  comparison 
to  the  generation  of  AE  in  the  war.i!  sector. 


V.  CONCLUSIONS 


The  purpose  of  this  study  is  to  determine  the  reasons 
for  the  growth  and  maintenance  of  Pacific  polar  lows.  In 
particular,  the  major  hypothesis  to  be  tested  is  that  the 
Pacific  polar  low  forms  through  baroclinic  instability.  To 
accomplish  this,  a  3-dimensional,  primitive  equation  model, 
including  convective  and  planetary  boundary  layer  parameter- 
izations  is  used  with  two  sets  of  initial  data.  The  first 
set  is  a  jet  profile  typical  of  mean  winter  conditions  and 
the  second  is  a  jet  profile  typical  of  conditions  present 
during  the  genesis  period  of  Pacific  polar  lows  (Mullen, 
1979).  A  linearized  version  of  the  model  is  used  for 
adiabatic  and  diabatic  simulations  of  the  Pacific  polar  low. 

The  linear  studies  indicate  that  medium-scale 
(wavenumber  12)  waves  are  most  unstable  for  mean  wi  iter 
conditions.  However,  wave-zonal  flow  interaction  and  surface 
friction  are  found  to  preferentially  damp  the  shortest  waves 
and  thereby  shift  the  most  unstable  mode  to  wavenumber  9. 
Energy  budgets  verify  that  this  wave  forms  through  barocl¬ 
inic  instability.  These  results  are  similar  to  those  found 
by  Gall  (1976a),  who  computed  the  growth  rate  spectrum  for 
wavenumbers  less  than  15  for  mean  winter  conditions. 

Although  this  study  includes  wavenumbers  15-30,  growth  rates 
of  these  waves  were  found  to  be  much  less  than  for  the 
longer  waves. 


The  genesis  jet  exhibits  a  peak  in  the  growth  rate 
spectrum  at  wavenumber  8,  which  also  grows  primarily  through 
baroclinic  instability.  However,  the  largest  growth  rate  is 
found  at  wavenumber  40.  This  disturbance  is  confined  to  the 
lowest  300  mb  of  the  troposphere  and  exhibits  a  helical 
circulation.  An  energy  budget  verifies  that  this  wave 
grows  through  baroclinic  instability  and  Kel vi n-Helmhol tz 
(vertical  shear)  instability.  Stone  (1966)  and  Orlanski 
(1968)  showed  that  Kel v i n-Hel mhol tz  instability  can  be 
significant  for  Richardson  numbers  less  than  1,  a  condition 
initially  satisfied  near  the  surface  in  the  genesis  jet 
case . 

Wavenumber  40  is  found  to  be  sensitive  to  the  initial 
conditions  and  physical  processes.  This  wave  exhibits  a 
relatively  short  period  of  linear  growth,  approximately  one 
day,  as  opposed  to  several  days  for  longer  waves.  It  is 
suggested  that  this  is  due  to  the  stabilization  of  the  flow, 
relative  to  wavenumber  40,  by  either  wavenumber  40  or  a 
combination  of  waves.  For  this  reason,  a  modified  genesis 
jet  is  created  by  allowing  pe rt urba t i ons  to  interact  with 
the  zonal  flow  for  84  h,  and  then  computing  the  new  zonal 
flow  by  averaging  all  84  h  forecast  fields  over  latitudinal 
circles.  In  this  modified  genesis  jet,  the  Richardson 
number  in  the  low-levels  near  the  jet  axis  ranges  from  2-5 
and  the  growth  rate  of  wavenumber  40  is  0.82  day-1.  This 
is  in  contrast  to  the  original  genesis  jet  which  exhibits 


Richardson  numbers  less  than  1  in  the  low-levels  underneath 
the  jet  axis  and  a  growth  rate  for  wavenumber  40  of  1.82 
day-1.  In  addition,  surface  friction  is  found  to  result  in 
a  further  reduction  in  the  growth  rate  of  wavenumber  40. 
However,  nonlinear  simulations  of  a  longer  wave  (wavenumber 
20)  indicates  the  presence  of  Kel  vi  n-Helinhol  tz  instability 
by  the  formation  of  bands  similar  to  that  found  in  wave- 
number  40  in  the  linear  experiments.  These  bands  are 
located  along  the  front  and  are  correlated  with  Richardson 
numbers  less  than  1  in  the  low  levels.  These  low  Richardson 
numbers  are  a  result  of  a  large  surface  sensible  heat  flux 
forcing  near-dry  adiabatic  lapse  rates  in  the  near-surface 
layers.  This  suggests  that  although  wave-zonal  flow  inter¬ 
action  and  surface  friction  act  to  damp  severely  the  growth 
of  the  barocl i ni c/Kel vi n-Helmhol tz  wave,  physical  processes 
can  create  favorable  conditions  for  its  existence. 

Growth  rates  for  the  modified  genesis  jet  are  signifi¬ 
cantly  different  from  the  original  genesis  jet  with  one 
growth  rate  maximum  at  wavenumber  8  and  another  of  equal 
size  at  wavenumber  20.  The  latter  disturbance  resembles  the 
Pacific  polar  low  described  by  Reed  (1979)  and  Mullen  (1979) 
in  that  it  forms  about  1000  km  north  of  the  jet  axis  with  a 
wavelength  of  1400  km.  However,  wavenumber  20  has  maximum 
amplitude  at  the  surface  and  extends  only  through  about  the 
lowest  300  mb  of  the  troposphere.  This  is  inconsistent  with 
observed  mature  Pacific  polar  lows  which  can  exhibit  height 


perturbations  at  300  mb  as  large  as  at  the  surface  (Reed, 
1979).  Energy  budgets  verify  that  wavenumber  20  forms 
through  baroclinic  instability.  The  magnitude  of  the  growth 
rates  found  in  this  study  for  wavenumbers  20  and  40  demon¬ 
strate  that  systems  with  horizontal  scales  of  the  Pacific 
polar  low  (and  smaller)  are  unstable  without  the  inclusion 
of  diabatic  processes. 

The  results  of  these  linear  studies  using  the  genesis 
jet  differ  from  those  of  Sardie  and  Warner  (1983)  (hereafter 
SW),  who  showed  that  the  maximum  growth  rate  for  the  genesis 
jet  is  1.29  day'1  at  wavenumber  10.  SW  found  no  evidence  of 
the  larger  growth  rates  for  shorter  wavelengths  as  found  in 
this  study.  Two  reasons  are  cited  for  this.  First,  SW  used 
a  quas i -geostrophi c  model  which  neglects  the  non-geostrophi c 
term  (u’w1  au/3p)  that  was  found  to  be  significant  in  the 
development  of  wavenumber  40.  Second,  the  model  used  by  SW 
contained  only  three  layers,  with  the  lowest  layer  being  200 
mb  thick,  whereas  the  model  used  here  contains  three  layers 
within  the  lowest  200  mb.  In  this  study,  for  linear,  adia¬ 
batic,  inviscid  conditions,  wavenumbers  20  and  40  are 
confined  to  the  lowest  250  mb,  which  suggests  that  the 
vertical  resolution  used  by  SW  is  insufficient  to  define 
these  features. 

Adiabatic  and  diabatic  nonlinear  simulations  are 
performed  for  wavenumber  20  of  the  modified  genesis  jet. 
Although  the  adiabatic  system  intensifies,  it  remains  a 


shallow  disturbance  as  in  the  linear  experiments  and  the 
deepening  rate  is  much  smaller  than  typically  observed  for 
Pacific  polar  lows.  The  inclusion  of  diabatic  effects 
dramatically  increases  both  the  disturbance  depth  and 
deepening  rate.  This  is  accomplished  by  enhancing  the 
temperature  perturbation  within  the  disturbance  such  that 
large  amounts  of  AE  are  generated  by  eddy  heat  transport. 
The  temperature  perturbation  is  increased  in  two  ways.  The 
first  is  through  latent  heat  release,  primarily  by  convec¬ 
tion,  in  the  warm  sector  of  the  storm.  Second,  large 
initial  air-sea  temperature  differences  result  in  a  strong 
surface  heat  flux  which  acts  to  maintain  the  low-level 
horizontal  temperature  gradient.  In  addition,  the  sensible 
heat  added  within  the  PBL  is  rapidly  mixed  to  mid-levels 
through  dry  convective  adjustment.  Rapid  intensification 
results  from  the  baroclinic  conversion  of  this  AE  to  KE. 

The  amount  of  AE  directly  generated  by  latent  heat  release 
in  the  warm  sector,  where  eddy  efficiency  factors  are 
positive,  is  about  an  order  of  magnitude  less  than  that 
produced  by  eddy  heat  transport.  Overall,  the  inclusion  of 
diabatic  effects  increased  the  kinetic  energy  of  the  polar 
low  by  over  an  order  of  magnitude. 

Significant  KE  is  generated  in  the  diabatic  simulation 
of  the  polar  low  by  Kel v i n-Hel mhol t z  instability.  Maximum 
values  of  this  instability  are  found  in  the  trough  asso¬ 
ciated  with  the  polar  low.  While  this  instability  does  not 


appear  to  be  significant  to  the  overall  development  of  the 
polar  low,  it  may  be  significant  in  the  formation  of  meso- 
scale  features  near  the  front  as  suggested  by  Bennetts  and 
Hoskins  (  1979  )  and  observed  by  Locatelli  e_t  aj_.  (  1982  ). 

The  simulated  polar  low,  with  diabatic  effects  included, 
is  found  to  resemble  a  simulated  synoptic-scale  extra- 
tropical  cyclone  (wavenumber  8),  as  observed  by  Locatelli  et 
al .  (1982).  Similarities  exist  in  patterns  of  wind,  tempera 
ture,  precipitation  and  surface  pressure.  The  effects  of 
diabatic  processes  on  wavenumber  8  is  found  to  be  similar  to 
that  for  the  polar  low. 

The  results  of  this  study  confirm  the  major  hypothesis 
that  the  Pacific  polar  low  is  a  baroclinic  disturbance. 
However,  diabatic  effects,  specifically  latent  heat  release 
in  the  warm  sector,  and  the  presence  of  an  initial  air-sea 
temperature  difference  of  at  least  5  K,  are  necessary  for 
the  polar  low  to  achieve  observed  deepening  rates.  This  is 
consistent  with  the  speculation  of  Sardie  and  Warner  (1983) 
that  moist  baroclinity  is  responsible  for  the  growth  of  the 
Pacific  pol ar  low. 

Research  should  be  continued  to  further  our  under¬ 
standing  and  prediction  capabilities  of  the  polar  low. 
Specifically,  our  understanding  of  convective  and  planetary 
boundary  layer  physics  must  continue  to  improve  since  they 
are  so  crucial  to  the  development  of  the  polar  low.  This 
can  be  accomplished  through  observational  studies.  More 


studies  are  needed  to  determine  the  precise  role  of  Kelvin- 
Helinholtz  instability  near  the  cold  front  and  its  possible 
role  in  the  formation  of  rainbands.  Also,  this  study  does 
not  address  either  the  Atlantic  polar  low  or  the  continental 
polar  low.  Studies  similar  to  this  one,  using  the  appro¬ 
priate  initial  conditions,  are  necessary  to  obtain  a  better 
understand i ng  of  these  systems. 


APPENDIX  A:  THE  FORECAST  MODEL 


A.  DYNAMICS 

The  forecast  model  uses  the  primitive  equations  in  the  a 
coordinate  system  after  Phillips  (1957), 
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The  vertical  structure  of  the  grid  is  shown  in  Fig.  A-l. 
All  variables  are  defined  at  the  middle  of  each  layer, 
except  tto ,  which  is  carried  at  the  interface  between  the 
layers.  The  sigma  levels  are  defined  in  Table  A-l. 

Table  A-l.  Model  levels. 

Level  Number  a 


1 

.05 

2 

.15 

3 

.24 

4 

.32 

5 

.40 

6 

.48 

7 

.57 

8 

.67 

9 

.76 

10 

.84 

11 

.92 

12 

.98 

The  horizontal  structure  of  the  grid  is  shown  in  Fig. 
A-2.  This  is  defined  as  scheme  C  by  Arakawa  and  Lamb 
(1977).  Here,  the  u-component  is  defined  one-half  grid 
interval  in  the  x-direction  from  the  mass  points  and 
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the  v-component  is  defined  one-half  grid  interval  in  the 
y-direction  from  the  mass  points.  The  grid  is  periodic 
in  x  and  has  free-slip  walls  on  the  y  boundaries. 

The  model  uses  the  split-explicit  time  integration 
scheme  developed  by  Madala  (1981).  First,  the  equations  are 
integrated  over  a  time  step  governed  by  the  meteorological 
and  slowest  gravity  modes.  The  solution  is  then  adjusted 
for  the  linear  terms  which  govern  the  fastest  gravity  modes. 
The  momentum  and  thermodynamic  equations  are  solved  with 
centered  time  and  space  differencing,  while  the  moisture 
equation  differencing  is  forward  in  time  and  upstream  in 
space.  Fourth  order  advection  is  used  for  the  momentum 
and  temperature.  Second  order  advection  is  used  for  the 
moi store . 

Fourth  order  diffusion  is  used  in  all  the  prognostic 
equations  except  for  *.  Second  order  diffusion  is  used  on 
the  first  interior  row.  This  diffusion  is  done  on  the  sigma 
levels,  which  can  lead  to  spurious  sources  and  sinks  of 
energy  when  applied  to  the  mass  and  moisture  fields.  To 
overcome  this  problem,  the  diffusion  is  applied  to  devia¬ 
tions  from  the  standard  atmosphere  for  the  temperature  and 
moisture.  The  diffusio.,  coefficient  is  computed  as: 

^  5  A  x / 3  on  1st  interior  row 

Ku  =  si. 6  Ax  on  2nd-4th  interior  rows  (A-9) 

) 

^1.25ax  on  all  other  interior  points 
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where  ax  is  the  grid  spacing  in  meters.  High  frequency  time 
oscillations  are  controlled  by  applying  the  Robert  (1966) 
time  filter  with  the  smoothing  coefficient  set  to  0.15. 

B.  PHYSICS 

1 .  Planetary  Boundary  Layer 

The  surface  fluxes  are  calculated  in  the  manner 
described  by  Deardorff  (1972).  First,  the  stability  of  the 
PBL  is  calculated  using  the  bulk  Richadson  number, 
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where  9V  is  the  vertical  potential  temperature,  subscript 
s  denotes  surface  values  and  subscript  m  denotes  mean  PBL 
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APPENDIX  C:  NONLINEAR  SIMULATIONS  OF  WAVENUMBER  8 


The  purpose  of  this  appendix  is  to  show  the  results  of 
nonlinear,  adiabatic  and  diaoatic  simulations  of  wavenumber 
3  for  the  genesis  jet  to  be  compared  to  the  results  for 
similar  simulations  of  wavenumber  20  presented  in  Chapter 
IV.  The  minimum  central  pressures  for  the  adiabatic  (G8A) 
and  diabatic  (G80)  experiments  are  given  in  Fig.  C-l.  The 
central  pressure  in  G8A  falls  from  1003  mb  at  72  h  to  968  mb 
at  120  h  with  rising  pressure  afterwards.  In  G8D ,  the 
central  pressure  falls  from  999  mb  at  72  h  to  944  mb  at 
138  h.  The  deepening  rate  is  similar  for  both  systems  from 
72-108  h.  However,  after  this  time,  G8D  exhibits  a  rapid 
decrease  in  central  pressure  while  the  central  pressure  in 
G8A  starts  to  rise  slowly.  The  reasons  for  this  pressure 
tendency  will  now  be  examined. 

The  sea-level  pressure  fields  for  60,  84,  108  and  132  h 
for  G8A  and  G30  are  shown  in  Figs.  C-2  and  C-3,  respective¬ 
ly.  At  60  h,  both  forecasts  exhibit  an  open  trough  pattern 
with  the  pressure  2-3  mb  lower  in  G8D.  In  addition,  the 
trough  extends  further  to  trie  southeast  in  G8D  with  the 
formation  of  a  second  center  in  this  trough.  The  extension 
of  this  trough  is  due  to  convection  in  the  warm  sector  of 
the  disturbance.  By  84  h,  G8A  has  attained  a  closed  circu¬ 
lation  wi  tn  a  central  pressure  of  996  mb,  whereas  G8D  has 


of  mass;  D£  (DZ)  is  the  dissipation  of  KE  (KZ)  by  surface 
friction;  GE  (GZ)  is  the  generation  of  AE  (AZ)  by  diabatic 
processes;  PAE  (PAZ)  is  the  change  in  AE  (AZ)  due  to  changes 
in  the  surface  pressure  and  PKE  (PKZ)  is  the  change  in  KE 
due  to  changes  in  the  surface  pressure. 
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(  B - 20  ) 


(  B  -  2 1 ) 


( B- 22  ) 


(  B  -  2  3  ) 


(  B  -  2  4  ) 


(  B  -  2  5  ) 


(  B  -  2  6  ) 


The  following  is  a  description  of  these  terms:  CA  is  the 
conversion  of  energy  between  AZ  and  A E ;  CE  is  the  conversion 
of  energy  between  AE  and  KE  by  vertical  r ed i s t r i b u t i on  of 
mass,  which  is  often  referred  to  as  barociinic  development; 
CK  is  the  conversion  of  energy  between  KZ  and  KE,  which  is 
often  referred  to  as  the  barotropic  term;  CZ  is  the  conver¬ 
sion  of  energy  between  KZ  and  AZ  by  vertical  red i s t r i b u t i o n 
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where  E  is  the  efficiency  factor,  defined  as, 


E  = 


p  k  _  p  k 

r 


( 8-6  ) 


The  reference  pressure,  P  ,  is  the  areal  mean  pressure  on  an 
isentropic  surface.  The  zonal  and  eddy  components  of  the 
kinetic  energy  are 
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respectively.  The  zonal  and  eddy  components  of  the 
available  potential  energy  are, 
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as  are  the  dynamics.  Therefore,  the  physical  parameteriza 
tions  only  need  to  be  calculated  at  some  multiple  of  the 
dynamics  time  step.  In  this  study,  the  non-convecti ve 
heating  and  the  cumulus  parameterization  are  applied  with 
a  2400s  time  step.  The  heating  and  moistening  rates  are 
distributed  uniformly  over  all  time  steps  until  these 
routines  are  called  again.  The  PBL  processes  and  the  dry 
convective  adjustment  are  called  every  dynamic  step. 


Moisture  is  condensed  until  saturation  is  reached.  All 
condensed  moisture  falls  into  the  next  layer  which  is  then 
checked  for  supersaturati on.  Precipitation  occurs  only  when 
the  air  is  saturated  from  the  cloud  to  the  ground. 

3 .  Dry  Convective  Adjustment 

Dry  convective  adjustment  is  applied  whenever  the 
lapse  rate  exceeds  the  dry  adiabatic  lapse  rate.  Tempera¬ 
tures  in  the  unstable  layer  are  modified  such  that  the  total 
potential  energy  is  conserved.  All  unstable  layers  are 
adjusted  subject  to  the  conditions 
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( A-36  ) 


and 
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(AT)  da  =  0  , 


( A- 37  ) 


where  aT  is  the  temperature  change  required  to  reduce  the 
lapse  rate  to  neutral.  These  equations  are  solved  as  a 
system  of  linear  equations. 

4 .  Incorporation  of  Physical  Processes 

Physical  processes  tend  to  evolve  relatively  slowly 
and  are  not  restricted  as  severely  by  time  step  limitations 


t 


» 


» .. 

9 


9 


142 


( A  -  3  3  ) 


sc 


=  IT  (o 


LCL 


aPBL )/2At 


where  At  is  the  time  step  for  the  physics  routine.  This 
acts  to  drive  the  PBL  halfway  down  to  the  LCL  in  the  absence 
of  other  effects. 

Cloud  temperature  (T  )  and  mixing  ratios  (qc)  are 
assumed  to  be  those  of  the  moist  adiabat  which  passes 
through  the  LCL.  The  partitioning  of  the  heating  and 
moistening  follows  Kuo  (1965).  A  value  b  is  defined  so 
that  a  fraction  bMt  is  used  to  moisten  the  environment  and 
(l-b)Mt  condenses  and  falls  instantaneously  as  rain.  The 
value  of  b  is  defined  as 

aT 

b  =  (aJ-cT)  /  (1'IRH)da  (A-34) 

B  T  aB 


where  IRH  is  the  integrated  relative  humidity  (0  £  IRH  _<  1) 
over  the  depth  of  the  cloud, 


IRH 


Mc 


( A- 3  5  ) 


Non-con vect i v e  precipitation  is  assumed  to  occur 
when  supersaturation  conditions  occur  at  a  grid  point. 
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specified  as  10°  at  all  times.  The  time  scale,  S,  in  ( A- 2  7 ) 
is  given  by 


o.75(eh-e0) 

deQ/dt 


( A- 30  ) 


when  (9h  -  80)  >  1°K,  and  by  Khakimov's  (1976)  equation 


0.67 

f 


(  A-  3 1 ) 


when  (0,  -  0  )  <  1  ° K . 

2 •  Preci pi tati on 

A  modified  version  of  the  Kuo  (1965)  cumulus 
parameterization  is  used.  This  modified  version  links  the 
convection  to  the  PBL.  First,  the  moisture  convergence,  Mt> 
is  cal cul ated  as 

Mt  -  l  v*^m7rVm)(1-aPBL)  +  psi^T)s  (A"32) 

Here,  Mt  the  moisture  convergence  occurring  solely  in  the 
PBL.  Convection  can  occur  when  Mt  >  4.0  X  10"®  gm  m-2s-l 
and  when  the  equivalent  potential  temperature  decreases  with 
height  from  the  PBL  to  the  first  model  layer  above  the  PBL. 
The  final  constraint  is  that  convection  cannot  occur  if  the 
LCL  is  above  the  PBL.  In  addition,  if  convection  is  occur¬ 
ring,  the  c 1 oud- i nd uc ed  subsidence  term  (term  C)  in  Eq. 
(A-23)  is  computed  as 
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The  entrainment  velocity  is  found  by  iterating  (A-24)  and 
(A-26).  The  first  term  on  the  right-hand  side  of 
(A-24)  describes  the  buoyant  production  of  turbulent 

kinetic  energy,  while  the  second  term  describes  the  mechan¬ 
ical  production  of  energy.  Stull's  equation  for  Wg  also 
contained  a  term  for  energy  loss  due  to  internal  gravity 
waves.  This  term  is  neglected  here.  Term  C  in  (A-23) 
accounts  for  cl oud- i nduced  subsidence.  The  computation  of 
this  term  is  described  later  in  this  Appendix. 

The  prognostic  equation  for  the  PBL  depth  in  stable 
conditions  follows  the  development  of  Niewstadt  and 
Tennekes  (1981). 


The  equilibrium  height  is  given  by 
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over  land  when  | d e  / d 1 1  >  0.1°K  hr-l  and 


h_  = 
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( A-29 


over  land  or  water  when  |dg  /dtj  <_  0.1°K  hr"l  after  Brost 


and  Wyngaard  (1978).  The  cross-isobar  angle  is  arbitrarily 


where  a  can  be  u,  w,  e  or  q.  By  making  these  assumptions, 
the  values  found  in  the  lowest  layer  of  the  model  are  always 
the  mean  quantities  in  the  PBL. 


The  height  of  the  PBL  in  unstable  conditions  is 
modeled  after  Stull  (1976).  The  prognostic  equation  for 
the  PBL  depth  (&P)  is 


irUP)  =  -  V-(WAp)  +  W  +  VJ 
oo  esc 

(a)  (b)  (c) 

Term  (a)  defines  the  effect  of  large-scale  motion  on  the 
PBL.  Term  (b)  is  the  entrainment  velocity  given  by 


( A- 23 ) 
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where 


Au 


( A- 2 5  ) 


A0vm  is  the  inversion  strength  at  the  PBL,  and 
is  the  depth  of  the  layer  of  negative  heat  flux  at  the  top 
of  the  PBL  , 


dl  = 


W  A0  h 
e  vm 


(w1 0 1  )  +  W  A0 

'  M  ’  ft 


( A- 2  6  ) 


respectively.  Once  C  and  C  are  known,  the  surface  stress 

u  o 


is  calculated  from 

2  2  2 
pu*  =  pum  Cu 

and  the  virtual  temperature  flux  is 


( A- 20  ) 


(-w'e;p  =  ceu*  ^  9  vm  -  evs)  *  (A-21) 

Deardorff  (1972)  proposed  extrapolation  equations  to 
determine  the  mean  PBL  quantities.  This  can  lead  to  serious 
problems  if  the  PBL  does  not  extend  at  least  through  the 
bottom  model  layer.  To  avoid  extrapolations  when  the  PBL 
is  shallow,  the  PBL  is  constrained  to  always  be  at  least  as 
deep  as  the  lowest  model  layer  (about  40  mb).  As  long  as 
the  lowest  layer  is  reasonably  shallow,  any  errors  intro¬ 
duced  should  be  relatively  small  since  such  a  shallow  PBL 
will  be  presumably  stable  and  the  surface  fluxes  should  be 
relatively  small. 

The  PBL  is  assumed  to  always  be  well-mixed.  Each 
time  the  model  physics  routines  are  called  (this  can  be  any 
integer  multiple  of  the  dynamic  time  step),  u,  v,  0  and  q  at 
all  model  levels  within  the  PBL  are  replaced  by  average 
values  throughout  the  PBL.  That  is, 
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and 


Ri c  =  3.05  ( A- 1 5  ) 

Rig  was  constrained  to  be  no  larger  than  90%  of  R i  ^ .  The 
surface  roughness,  1Q  is  a  specified  constant  over  land  and 

ZQ  =  0.01  u*2  (  A- 1 6  ) 

over  water.  The  value  of  ZQ  over  water  is  constrained  to  be 
no  less  than  2.43  5  X  1  0“4  m  (suggested  by  Dr.  S.  Burk, 
NEPRF). 

For  unstable  conditions  (Ri_  <  0),  the  friction 

D 

coefficient  is 

Cu  =  [Cun_1  -  25  exP  (°-26C  -  0 . 03  C2 ) ] " 1  ,  (A-17) 

and  the  heat  transfer  coefficient  is 

C„  =  (Cfl  _1  +  C  _1  +  C  _1 )-1  ,  ( A- 1 8 ) 

9  9n  u  un  '  v 

where 

C  =  log1Q  (-RiB)  -  3.5  ( A - 1 9 ) 

and  C.._  and  Can  are  given  by  (  A  - 1  3  )  and  (  A  -  1 4  )  , 
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Fig.  C-l.  Minimum  central  pressure  for  G8A  (solid)  and  G8D  (dashed). 
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Fig.  C- 3 . 


Similar  to  Fig.  C-2  except  for  G8D. 


deepened  to  991  mb.  By  108  h,  both  systems  have  intensified 
further  with  G8A  near  the  mature  stage.  At  132  h,  G8A 
remains  at  nearly  the  same  central  pressure  while  G8D  has 
continued  to  deepen  and  has  reached  the  mature  stage  with  a 
central  pressure  of  947  mb. 

Cross-sections  of  height  perturbations  through  G8A  and 
G8D  at  84  h  are  given  in  Fig.  C-4.  Both  systems  exhibit 
maximum  perturbations  at  the  surface  and  in  the  upper  levels 
with  a  westward  tilt  with  height.  The  upper  level  height 
perturbation  is  stronger  in  G80  than  in  G8A,  which  shows 
that  diabatic  effects  tend  to  strengthen  the  upper  level 
circulation.  This  is  also  found  to  be  the  case  for  the 
polar  low  (cf,  Fig.  4-5).  Unlike  G20A,  which  is  strongly 
damped  in  the  vertical,  the  upper  level  perturbation  in  G8 A 
is  more  than  half  the  magnitude  of  the  surface  perturbation. 
This  again  illustrates  that  shorter  waves  tend  to  be  shallow 
disturbances  as  found  in  Chapter  III. 

Cross-sections  of  vertical  motion,  temperature  and  the 
trough  axis  for  G8A  and  G8D  at  84  h  are  shown  in  Fig.  C-5. 
The  adiabatic  case  exhibits  a  warm  perturbation  extending 
from  the  surface  to  400  mb.  The  center  of  the  warmest  air 
tilts  eastward  with  height  from  the  surface  to  500  mb,  and 
westward  above  500  mb.  A  center  of  upward  vertical  motion 
occurs  just  east  of  the  surface  low  in  the  warmest  air, 
which  implies  an  increase  of  KE  over  a  deep  layer  due  to 
this  thermally  direct  circulation.  This  is  in  contrast  to 


the  polar  low,  G20A  (Fig.  4-2),  which  exhibits  very  little 
tilt  to  the  warm  air  perturbation  from  the  surface  to  750  mb 
with  a  strong  westward  tilt  above  750  mb. 

Significant  differences  exist  between  the  cross-sections 
of  G8A  and  G8D.  First,  the  center  of  upward  vertical  motion 
and  the  warmest  air  in  G80  are  found  about  750  km  east  of 
their  respective  positions  in  G8A.  This  is  a  result  of  the 
convection  taking  place  ahead  of  the  surface  low  as  noted 
above.  Second,  in  G8D,  a  narrow  band  of  upward  motion 
within  a  weak  warm  air  perturbation  is  found  in  the  cold  air 
about  1000  km  behind  the  position  of  the  surface  low.  The 
surface  reflection  of  this  feature  is  the  weak  pressure 
trough  in  the  northwest  flow  in  G80  at  84  h  (Fig.  C-3b). 

This  pressure  trough  is  a  weak  polar  low  and  its  formation 
is  related  to  an  interaction  between  convection  and  the  PBL. 
Forecast  PBL  depths  in  the  cold  air  behind  the  trough  are 
relatively  large  because  the  maximum  air-sea  temperature 
differences  are  in  this  area.  Some  convection  (about  5  mm/ 

12  h)  is  forecast  behind  the  surface  low  in  the  cold  air. 
This  results  in  a  decrease  in  the  PBL  height  at  the  convec¬ 
tive  points  due  to  cl oud- i nduced  subsidence  (Eq.  A-33).  At 
the  points  where  the  PBL  is  deep,  relatively  strong  PBL 
winds  exist  due  to  the  mixing  of  the  winds  throughout  the 
PBL.  At  the  convective  points,  where  the  PBL  is  shallow, 
relatively  light  winds  exist  since  the  stronger  winds  aloft 
are  not  being  mixed  into  the  PBL.  This  sets  up  a  zone  of 


convergence  in  the  PBl  on  the  western  edge  of  convective 
areas  which  further  enhances  convective  forcing  of  this 
system.  While  this  process  is  capable  of  producing  small- 
scale,  cold  air  disturbances  in  this  model,  it  is  felt  that 
this  is  more  of  a  numerically  induced  feature  as  opposed  to 
a  physically  realistic  situation.  Certainly,  not  all 
convection  drives  the  PBL  halfway  down  to  the  LCL  every  half 
hour  as  is  assumed  here.  However,  large  variations  in  PBL 
depth  over  a  small  area  in  the  pretence  of  strong  vertical 
wind  shear  should  be  considered  as  a  possible  mechanism  to 
trigger  polar  low  development. 

The  total,  convective  and  stable  precipitation  accumu¬ 
lated  in  G8D  during  the  12  h  period  ending  at  108  h  are 
shown  in  Fig.  C-6.  About  26  mm  of  precipitation  is  forecast 
to  the  north  of  the  disturbance  with  a  maximum  of  38  mm  to 
the  southeast  of  the  low  along  the  cold  front.  Only  about 
40%  of  the  northern  maximum  and  2/3  of  the  southern  maximum 
in  precipitation  is  convective.  Therefore,  the  percentage 
of  convective  precipitation  is  smaller  in  G8D  than  in  the 
polar  low,  G20D. 

The  effect  of  the  diabatic  processes  can  be  seen  in 
Figs.  C-7  and  C-8  in  which  the  850,  700  and  500  mb  height 
and  temperature  fields  are  given  for  G8A  and  G80,  respec¬ 
tively.  As  with  the  polar  low  in  Chapter  IV,  a  noticeable 
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Fields  of  height  and  temperature  at  (a)  850  mb 
(b)  700  mb  and  (c)  500  mb  at  108  h  for  G8 A . 
Contour  interval  for  height  is  60  m  in  (a)  and 
(b)  and  120  m  in  (c).  Grid  placement  as  in 


Fig.  C-8.  Similar  to  Fig.  C-7  except  for  G8D. 


warming  occurs  ahead  of  the  developing  wave  due  to  the 
release  of  latent  heat  in  this  area,  which  reduces  the 
effect  of  adiabatic  cooling  through  ascent. 

The  surface  sensible  and  latent  heat  fluxes  associated 
with  G8D  at  108  h  are  shown  in  Fig.  C-9.  Maximum  values  of 
the  fluxes  are  found  behind  the  developing  wave  in  the 
colder  air,  as  noted  with  the  polar  low  (cf,  Fig.  4-11). 
However,  unlike  the  polar  low,  downward  (negative)  fluxes 
occur  ahead  of  G80  in  the  warm  sector.  This  results  from  a 
large  low-level  temperature  perturbation  (Fig.  C-8a)  in  G8D 
which  reverses  the  sign  of  the  air-sea  temperature  differ¬ 
ence  in  the  warm  sector. 

KE  budgets  for  G8A  and  G8D,  integrated  over  the  forecast 
domain,  are  presented  in  6  h  intervals  in  Tables  C-l  and 
C-2,  respectively.  The  maximum  KE  attained  by  G8 A ,  6.7  x 
1 0 5  j  m-2>  i$  about  40%  less  than  the  maximum  attained  by 
G8D,  11.1  x  10^  J  m-2.  This  is  in  sharp  contrast  to  the 
results  in  Chapter  IV  that  the  diabatic  effects  led  to  an 
order  of  magnitude  increase  in  KE  for  G20D.  The  reason  why 
G8D  does  not  exhibit  a  similar  increase  in  magnitude  over 
G8A  can  be  seen  from  Fig.  C-10,  which  gives  the  time-height 
cross-sections  of  KE  for  G8A  and  G80.  For  the  polar  low 
discussed  in  Chapter  IV,  the  adiabatic  run  exhibits  a  maxi¬ 
mum  of  KE  near  the  surface  with  negligible  amounts  above 
600  mb.  However,  G8A  exhibits  a  maximum  KE  of  5  x  1  0  4  J  m-2 
near  350  mb  with  a  secondary  maximum  of  3  x  104  J  m-2  near 


Table  C-l.  KE  budgets  for  G8A.  KE  in  units  of  10^  Jm"^ 
All  other  budget  terms  in  units  of  Wm_2. 


Time  (hr) 

KE 

3  KE/a  t 

CE 

HCK 

VC  K 

DE 

60 

11.9 

.29 

.31 

-.06 

.11 

-.08 

66 

19.9 

.49 

.51 

-.11 

.18 

-.13 

72 

32.9 

.77 

.84 

-.19 

.29 

-.22 

78 

53.0 

1.16 

1.35 

-.34 

.44 

-.37 

84 

83.0 

1.68 

2.10 

-.57 

.62 

-.62 

90 

125.5 

2.33 

3.12 

-.82 

.81 

-1.03 

96 

183.6 

3.08 

4.36 

-1.03 

.99 

-1.63 

102 

258.4 

3.78 

5.73 

-1.26 

1.16 

-2.48 

108 

346.9 

4.22 

7.12 

-1.71 

1.37 

-3.60 

114 

440.8 

4.21 

8.16 

-2.33 

1.64 

-4.85 

1  20 

528.8 

3.68 

8.46 

-2.82 

1.85 

-5.91 

126 

600.0 

2.77 

7.98 

-3.09 

1.92 

-6.55 

132 

648.3 

1.67 

6.94 

-3.03 

1.78 

-6.72 

138 

672.1 

.48 

5.66 

-2.71 

1.54 

-6.71 

144 

669.1 

-.78 

3.93 

-2.23 

1.27 

-6.51 

150 

638.6 

-1.80 

2.22 

-1 . 94 

1.05 

-5.71 

in  Fig.  C - 1 6 .  Although  G3 A  suggests  a  maximum  near  450  mb 
and  another  near  800  mb,  both  forecasts  exhibit  large  values 
throughout  the  900-300  mb  layer.  A  similar  pattern  is  found 
for  G20D  except  over  a  slightly  shallower  layer. 

Noticeable  differences  exist  in  the  vertical  eddy  heat 
transport  term,  VCA,  between  G8A  and  G8D.  In  the  adiabatic 
run,  VCA  is  positive  during  the  developing  stage  up  to  102 
h,  which  indicates  that  VCA  is  a  source  of  AE.  After  102  h, 
VCA  is  negative  and  tends  to  damp  the  wave.  However,  values 
of  VCA  are  generally  4-10  times  less  than  the  largest  term, 
HCA.  In  G8D,  VCA  is  positive  through  90  h.  While  G8D  is 
deepening  rapidly  from  96-144  h,  VCA  is  negative  (except  at 
120  h)  and  is  very  small  in  magnitude.  Even  though  VCA  is 
negative  during  this  time,  the  values  are  not  as  large  in 
magnitude  as  in  G8A,  and  therefore  more  AE  is  available  in 
G8D.  Time-height  cross-sections  of  VCA  for  G8A  and  G8D  are 
given  in  Fig.  C - 1 7 .  Generally,  negative  values  are  found 
above  500  mb.  Below  500  mb,  G8A  is  dominated  by  positive 
values  prior  to  108  h  and  negative  values  after  108  h.  The 
same  tendency  exists  in  G8D  except  that  positive  values 
persist  at  all  times  below  800  mb. 

The  destruction  of  AE  by  VCA  in  G8A  and  G8D  during  the 
growth  stage  is  different  from  that  found  in  the  polar  low. 
In  G20D,  generation  of  AE  through  VCA  during  the  developing 
stage  was  found  to  be  approx i matel y  the  same  magnitude  as 
term  HCA.  The  reason  for  this  appears  to  be  the  area  in 
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Fig.  C-15.  Time-height  cross-section  of  GE  for  G8D. 
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and  G3D  occur  in  the  700-500  mb  layer  from  120-132  h  (Fig. 

C - 1 1 )  .  This  is  at  the  same  time  and  location  in  which  the 
large  increases  in  AE  are  observed,  which  suggests  that  the 
disturbance  uses  the  additional  AE  for  baroclinic  develop¬ 
ment  of  the  storm. 

The  generation  of  AE  by  diabatic  processes,  GE,  attains 
maximum  values  from  114-138  h.  This  is  consistent  with  the 
large  increases  in  AE  noted  during  this  time.  The  time- 
height  cross-section  of  GE  for  G8D  is  given  in  Fig.  C - 1 5 . 

The  largest  values  from  114-138  h  are  found  in  the  700- 
500  mb  layer.  This  is  attributed  to  the  latent  heat  release 
primarily  in  the  warm  sector  ahead  of  the  low  as  discussed 
above.  Negative  values  of  GE  in  the  lowest  100  mb  are  the 
result  of  positive  (negative)  sensible  heat  fluxes  occurring 
in  the  cold  (warm)  air  behind  (ahead  of)  the  low.  This 
pattern  of  GE  resembles  that  for  G20O  (Fig.  4-17). 

The  horizontal  eddy  heat  transport,  HCA,  is  positive  for 
both  G8A  and  G80,  which  indicates  that  it  is  a  source  of  AE. 
During  the  developing  stages  of  G8A  and  G8D,  HCA  is  the 
largest  term  in  the  AE  budget.  Values  of  HCA  are  up  to  two 
times  larger  in  G8D,  which  suggests  that  the  diabatic 
processes  have  altered  the  overall  temperature  distribution 
of  the  disturbance  so  as  to  increase  the  horizontal  eddy 
heat  transport.  This  was  shown  in  Chapter  IV  to  play  a 
significant  role  in  the  intensification  of  the  polar  low. 
Time-height  cross-sections  of  HCA  for  G8A  and  G8D  are  given 


Table  C-4.  Similar  to  Table  C-3  except  for  G8D 


Time  (hr) 

AE 

3  AE/3  t 

CE 

GE 

HCA 

VCA 

60 

11.7 

.37 

.70 

.10 

.53 

.52 

66 

22.2 

.63 

1.42 

.32 

1.05 

.89 

72 

38.8 

.88 

2.01 

.52 

1.95 

.76 

78 

60.1 

1.08 

2.69 

.33 

3.13 

.66 

84 

85.8 

1.20 

4.28 

1.09 

4.54 

.78 

90 

111.9 

1.36 

5.28 

1.38 

5.82 

.07 

96 

144.5 

1.70 

5.77 

1.48 

7.51 

-.36 

102 

185.5 

2.39 

7.39 

1.88 

9.65 

-.60 

108 

247.8 

3.13 

9.58 

2.51 

12.33 

-.44 

114 

320.9 

3.30 

12.14 

3.08 

14.48 

-.11 

120 

390.5 

2.47 

13.84 

3.39 

14.97 

.22 

126 

427.7 

.82 

13.23 

2.88 

13.95 

-.03 

132 

426.1 

-.92 

12.81 

3.05 

11.99 

-.90 

138 

387.8 

-2.05 

10.63 

2.97 

8.79 

-1.09 

144 

337.6 

-2.22 

7.11 

2.02 

5.18 

-.43 

150 

291.9 

-1.92 

3.57 

1.09 

1.68 

.37 
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Table  C-3.  AE  budgets  for  G8A .  AE  in  units  of  103 
All  other  budget  terms  in  units  of  Wm'2 


Time  (hr) 


AE  3AE / 3t  CE  HCA  VCA 


cross- sect i ons  of  VCK  for  G8A  and  G80  are  given  in  Fig. 

C - 1 3 .  Both  disturbances  exhibit  maximum  values  of  VCK  of 
approximately  the  same  magnitude  between  850-900  mb  and  near 
400  mb. 

Frictional  dissipation,  DE,  acts  as  a  sink  of  KE  for 
both  forecasts.  DE  is  the  largest  term  in  both  G8A  and  G8D 
once  the  disturbance  reaches  its. maximum  intensity.  This  is 
in  contrast  to  G20D,  which  exhibited  larger  values  of  HCK 
than  DE  during  the  decaying  stage  of  the  polar  low.  How¬ 
ever,  the  trends  here  are  similar  since  for  G8A,  DE  is  2.2 
times  larger  than  HCK  at  132  h  while  for  G8D ,  DE  is  1.7 
times  as  large  as  HCK  at  the  same  time.  Therefore, 
barotropic  damping  is  larger  relative  to  DE  when  diabatic 
effects  are  included. 

/TE  budgets  integrated  over  the  forecast  domain  are 
presented  in  6  h  intervals  for  G8A  and  G8D  in  Tables  C-3  and 
C-4,  respectively.  As  was  the  case  for  the  polar  low,  G8D 
exhibits  more  AE  than  G8A.  The  maximum  values  for  G8A,  2.2 
x  1 0  5  j  m-2,  and  for  G8D,  4.3  x  1 0  5  J  m-2,  occur  at  126  h. 
Time-height  cross-sections  of  AE  for  G8A  and  G8D  are  given 
in  Fig.  C-14.  Although  G8D  exhibits  larger  values  of  AE  at 
all  levels,  the  most  significant  increase  over  G8A  occurs  in 
the  mid-levels  from  700-500  mb.  AE  values  for  G8D  in  this 
layer  are  three  times  larger  from  120-132  h  when  the  central 
pressure  fell  7  mb  in  G8D  and  rose  2  mb  in  G8A.  It  was 
noted  earlier  that  the  largest  differences  in  CE  between  G8A 


(8 w)  3anss3yd 


in  G8A.  The  vertical  structure  of  CE  is  very  similar 
between  the  diabatic  simulations  of  G20  ( c f ,  Fig.  4-13) 
and  G8 . 

The  horizontal  barotropic  term,  HCK,  is  a  sink  of  KE  in 
both  forecasts.  Maximum  values  of  HCK  are  found  at  the 
mature  stages  of  G8A  and  G8D.  Time-height  cross-sections  of 
HCK  are  given  in  Fig.  C - 1 2 .  The  largest  values  of  HCK  are 
found  in  the  upper  levels  associated  with  the  jet.  Also  in 
G8D,  positive  values  of  HCK  are  found  below  600  mb  from  96- 
114  h.  The  reason  for  this  is  the  southeastward  extension 
of  the  developing  trough  forced  by  convection  as  seen  in  the 
surface  pressure  field  for  G8D  (Fig.  C-3).  This  pattern 
implies  a  positive  correlation  between  u'  and  v'  and  hence 
a  conversi on . of  KE  to  KZ  by  HCK. 

The  vertical  component  of  CK,  VCK,  is  positive  through¬ 
out  both  forecasts,  and  therefore  represents  a  source  of  KE. 
As  with  the  polar  low  forecasts,  VCK  is  a  maximum  6-12  h 
after  CE  reaches  its  maximum  value  near  the  mature  stage  of 
the  cyclone.  This  should  be  expected  since  u*  and  u'  are 
largest  then.  The  fact  that  VCK  attains  its  largest  values 
at  or  just  after  maximum  storm  intensity  suggests  that  this 
process  is  a  symptom  of  the  baroclinic  disturbance  rather 
than  a  cause,  as  is  also  the  case  for  HCK.  Therefore,  both 
G20  and  G8  are  considered  to  be  disturbances  which  form 
through  conventional  baroclinic  instability.  Time-height 
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100 


Contour  interval  is  20  x  10"^  W  m~ ^ .  Heavy  so 
zero  contour,  light  solid  line  is  positive  con 
dashed  line  is  negative  contour. 


900  mb.  Therefore,  the  reason  G20D  realizes  such  a  large 
increase  of  KE  when  diabatic  effects  are  included  is  that 
G20A  is  a  very  shallow  system  and  the  dom a i n- i nteg rated  KE 
is  quite  small.  With  diabatic  effects,  the  polar  low 
extends  through  a  large  part  of  the  troposphere  and  yields 
a  much  larger  domai n- i ntegrated  KE.  G8A  already  extends 
through  the  troposphere  and  the  only  other  way  to  increase 
the  domain-integrated  KE  is  to  increase  the  circulation  of 
the  system.  This  occurs  in  G8D  and  results  in  KE  maxima  of 
9  x  104  J  m-2  at  350  mb  and  6  x  104  j  m-2  at  900  mb. 

The  baroclinic  development  term,  CE,  is  the  largest  term 
in  the  budgets  for  both  G8A  and  G8D  through  132  h.  During 
the  developmental  stages  prior  to  120  h,  CE  is  1-1/2  to  2 
times  larger  in  G8D  and  is  the  largest  term  in  the  KE 
budgets.  Time-height  cross-sections  of  CE  are  given  in 
Fig.  C - 1 1 .  Both  forecasts  exhibit  maximum  conversion  of  AE 
to  KE  at  700  mb,  which  is  near  the  level  of  maximum  vertical 
motion.  Large  differences  exist  between  the  vertical  struc¬ 
ture  of  CE  for  G8A  and  G20A.  G20A  exhibits  maximum  CE  near 
900  mb  with  negative  values  above  700  mb.  This  is  a  result 
of  the  warm  air  perturbation  tilting  sharply  westward  above 
750  mb.  However,  the  warm  air  perturbation  in  G8A  tilts 
eastward  with  height  up  to  500  mb  during  the  development 
stage.  This  results  in  positive  values  of  CE  up  to  300  mb 


me-height  cross 
)  G80.  Contour 


Fig.  C-  1 6 


which  the  low  develops.  VCA  is  comprised  of  the  four 
terms:  a)  w'a'E/Cp;  b)  a/ Cp  ;  c)  (  wT'  8Ef/3p);  and  d) 

( -  oi'  E  1  3T/  3p ) .  Term  a  is  generally  an  order  of  magnitude 
less  than  the  other  terms,  so  VCA  is  a  balance  between  term 
b,  which  is  typically  less  than  zero,  and  terms  c  and  d, 
which  are  typically  positive.  The  relative  strength  of  each 
term  governs  the  sign  of  VCA.  In  G20D,  term  b  tends  to  be 
relatively  small  since  this  disturbance  forms  250-500  km 
farther  to  the  north  than  in  G8  and  therefore  a  is  smaller. 
This  reduces  the  magnitude  of  term  b  and  VCA  is  positive. 

The  results  of  the  G8  simulations  presented  here  illus¬ 
trate  the  following  points.  First,  diabatic  processes  serve 
to  intensify  an  ext  rat ro p i c a  1  ,  synopt i c- seal e  disturbance  in 
a  similar  manner  to  the  Pacific  polar  low.  Specifically, 
the  diabatic  processes  serve  to  alter  the  temperature  struc¬ 
ture  of  both  G8  and  G20  such  that  eddy  heat  transport,  and 
therefore  values  of  AE ,  are  increased.  The  direct  genera¬ 
tion  of  AE  by  diabatic  processes  (GE)  is  found  to  be  at 
least  5  times  smaller  than  the  eddy  heat  transport  in  both 
G8  and  G20.  Second,  patterns  of  surface  pressure,  precipi¬ 
tation,  wind  and  temperature  for  a  Pacific  polar  low  are 
found  to  resemble  those  of  synoptic-scale  oceanic 
di sturbances . 
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